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ABSTRACT 

We present the analysis of the full bi-dimensional optical spectral cube of the nearby spiral 
galaxy NGC 5668, observed with the PPAK IFU at the Calar Alto observatory 3.5m telescope. 
We make use of broad-band imaging to provide further constraints on the evolutionary history of 
the galaxy. This dataset will allow us to improve our understanding of the mechanisms that drive 
the evolution of disks. We investigated the properties of 62 H II regions and concentric rings 
in NGC 5668 and derived maps in ionized-gas attenuation and chemical (oxygen) abundances. 
We find that, while inwards of r~ 36" ~ 4.4kpc ~ 0.36 (^p-) the derived O/H ratio follows the 
radial gradient typical of spiral galaxies, the abundance gradient beyond r~ 36" flattens out. 
The analysis of the multi-wavelength surface brightness profiles of NGC 5668 is performed by 
fitting these profiles with those predicted by chemo-spectrophotometric evolutionary models of 
galaxy disks. From this, we infer a spin and circular velocity of A=0.053 and v c =167kms , 
respectively. The metallicity gradient and rotation curve predicted by this best-fitting galaxy 
model nicely match the values derived from the IFU observations, especially within r ~ 36". The 
same is true for the colors despite of some small offsets and a reddening in the bluest colors beyond 
that radius. On the other hand, deviations of some of these properties in the outer disk indicate 
that a secondary mechanism, possibly gas transfer induced by the presence of a young bar, must 
have played a role in shaping the recent chemical and star formation histories of NGC 5668. 



Subject headings: Galaxies: abundances — Galaxies: 
Galaxies: ISM — Galaxies: kinematics and dynamics 
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evolution — Galaxies: individual: NGC 5668 — 
-Techniques: spectroscopic 

1. Introduction 

The evolution of galactic disks is one of the 
most important and yet not fully understood top- 
ics in extragalactic astronomy. Despite signifi- 
cant progress in the recent past regarding our 
understanding of the history of both the thick 
and thin disks, important questions remain unan- 
swered: How old are the disks seen in the spiral 
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Fig. 1. — Left: SDSS optical primage (from 4458 to 7706 A) of NGC 5668 centered at RA(2000) 14 h 33 m 24.3 s and Dec(2000) 
+04°33'24.3". North is up and East is to the left. Plate scale is 0.396 arcsec/pix. Right: Synthetic false-color image obtained 
from the PPAK data cube and the response curves of the SDSS u, g and r filters. Note that the luminous spot in pointing 
7 is caused by cross-talk effects from a nearby (on the CCD) field star. Plate scale is 1 arcsec/pix in this case. Center and 
orientation are the same as in the left panel. 



galaxies today? How did they chemically evolve? 
Are they growing inside-out, as proposed to ex- 
plain the color and metallicity gradients of our 
own Milky Way? Do they have an edge? How 
efficient is stellar radial diffusion? 
Until recently, the study of the properties of spi- 
ral disks has been limited to broad-band imag- 
ing data and/or long-slit spectroscopy, which has 
severely limited the reach of previous works. Ra- 
dial color profiles have been widely used to probe 
the inside-out scenario of disk evolution (Bell & de 
Jong 2000; MacArthur et al. 2004; Pohlen & Tru- 
jillo 2006; Muhoz-Mateos et al. 2007). However, 
radial variations of metallicity and dust extinction 
conspire to create color gradients similar to those 
due to changes in stellar age alone. Analysis of 
color-magnitude diagrams allow for more robust 
measurements of stellar ages as a function of ra- 
dius, but these kind of studies are currently lim- 
ited to very nearby galaxies where individual stars 
can be resolved (de Jong et al. 2007; Barker et al. 
2010; Gogarten et al. 2010). Long-slit spectra 
overcomes the problems of broad-band photom- 
etry, but at the expense of focusing along a pre- 



defined spatial axis (Yoachim & Dalcanton 2008). 
Our effort is committed to add another dimension 
to the study of nearby spiral galaxies thanks to the 
use of wide-field integral-field spectrometers (IFS) 
combined with images covering a wide wavelength 
range from the ultraviolet to the infrared. This 
paper presents the analysis of the IFU data taken 
for NGC 5668 within a more ambitious project 
aimed to map a sample of a dozen nearby galax- 
ies included in the HGRS (Hcrschel Galaxy Ref- 
erence Survey; Boselli et al. (2010)) using mosaics 
obtained with two of the largest IFU available to 
date: PPAK (Pmas fiber PAcK) IFU of PMAS 
(Potsdam Multi- Aperture Spectrophotometer) at 
the CAHA 3.5-m and VIMOS (Visible MultiOb- 
ject Spectrograph) at ESO-VLT. 

The evolution of galaxy disks is a complex pro- 
cess as many are the mechanisms that might al- 
ter their photometric, chemical, and kinematical 
properties: minor mergers, formation of bars or 
rings, density waves, stellar diffusion, gas in-fall, 
and, in the case of galaxies in clusters, also ram 
pressure striping and galaxy harassment. Despite 



this rather discouraging scenario, a picture for the 
formation of disks has emerged in recent years 
that attempts to explain most of their observa- 
tional properties: the inside-out formation sce- 
nario. Although this scenario has been around 
for a long time now, it has been only recently 
that inside-out growth of galaxy disks has been 
applied and quantified beyond the MW. Both 
the extinction-corrected color gradients in nearby 
galaxies, Munoz-Mateos et al. (2007), and weak 
dependence of the mass-size relation with redshift 
(Trujillo et al. 2004, Barden et al. 2005, Trujillo 
et al. 2006) support an inside-out scenario for the 
evolution of disks. However, there are several ob- 
servational results that conflict or, at least, cannot 
be explained by this otherwise elegant scenario. In 
this regard, many nearby spiral disks show that, 
while the color and mctallicity gradients of the 
inner disk are compatible with the inside-out sce- 
nario, beyond the SF-threshold radius (Martin & 
Kennicutt 2001) the gradients flatten and even re- 
verse, by getting redder toward the outer regions 
(Bakos et al. 2008; Azzollini et al. 2008; Vlajic 
et al. 2009). In the specific case of the chemical 
abundances, the models based on the standard 
inside-out scenario of galaxy disk formation pre- 
dict a relatively quick self enrichment with oxygen 
abundances 1/10 the Solar value after only 1 Gyr 
of evolution (Boissier & Prantzos 1999). These 
models also predict an almost universal negative 
metallicity gradient once this is normalized to the 
galaxy optical size d[0/H] of ~ —0.8 dex/R 2 5; 
(Prantzos & Boissier 2000). Regarding this, the 
few observational abundance measurements ob- 
tained to date in the outer edges of disks show a 
flattening or even an increase in the abundance 
toward the outermost parts of the disk. This is 
also true in the case of the extended UV disks 
recently discovered by the GALEX satellite (Gil 
dc Paz et al. 2005; Gil dc Paz et al. 2007; Thilker 
et al. 2007), which show oxygen abundances that 
are rarely below 1/10 (Solar value). In this regard, 
Bresolin et al. (2009a) found that the normalized 
radial metallicity gradient in M83 changes from 
-0.030 dex kpc" 1 in the inner disk to -0.005 
dex kpc -1 in the outer parts, see also Bresolin 
et al. (2012). We also refer the reader to Yoachim 
et al. (2010) for another recent example in this 
same regard. Similar results were obtained on the 
mctallicity gradient of the outer disk of NGC 300 



from single-star CMD analysis (Vlajic et al. 2009). 

In this work we present the pilot study of the 
full bi-dimensional spectral cube of the nearby 
spiral galaxy NGC 5668, (Fig. 1), obtained with 
the PPAK Integral Field Unit at the Calar Alto 
(CAHA) observatory 3.5 m telescope. Despite the 
relatively modest collecting area of the CAHA 3.5- 
m telescope, the broad spectral coverage (from 
3700 to 7000 A) and adequate spectral resolution 
(R = 500) of PPAK with the V300 grating make 
this instrument one of the best tools for study- 
ing stellar populations, dust content, and physical 
conditions of the gas (temperature, density, chem- 
ical abundances) in spatially-resolved galaxies. 

2. The nearby spiral galaxy NGC 5668 

NGC 5668 is a nearly face-on late-type spiral 
galaxy classified as an Sc(s)II-III type according to 
Sandage & Tammann (1987) and as an SA(s)d by 
de Vaucouleurs et al. (1991). There is a weak bar 
or oval inner structure 12" in size visible on the op- 
tical image, which leads to a small shoulder in the 
surface brightness profile published by Schulman 
et al. (1994). The outer disk (beyond R = 100") 
is slightly asymmetric and more extended towards 
the North. For this work we adopted a distance 
of ~ 25 Mpc, (m — M) = 31.99 mag, assuming 
a cosmology with H = 73kms~ 1 Mpc~ 1 , S! ma ttcr 
= 0.27, ^vacuum = 0.73, for a recession velocity 
of 1813 kms -1 corrected to the reference frame 
defined by the 3K Microwave Background Radi- 
ation (Fixsen et al. 1996) 1 . A broad-band opti- 
cal image of this galaxy from the Sloan Digital 
Sky Survey (York et al. 2000) is shown in Fig. 1. 
The galaxy presents an inclination of 18° and a 
total magnitude B = 12.13 ± 0.03 mag (Schul- 
man et al. 1996). The total mass within the opti- 
cal disk was estimated to be 5.7xl0 10 M Q (Schul- 
man et al. 1996) 2 . NGC 5668 has been found to 
host a number of HVCs (High Velocity Clouds) 
by Schulman et al. (1996) using the Arecibo tele- 
scope. These authors detected high velocity wings 
in the line shape, which were attributed to HVCs 
in the galaxy. High velocity clouds and High- 



1 Source: NASA/IPAC Extragalactic Database 
(http: / /nedwww.ipac. caltech.edu/) 

2 They adopted a distance of 30 Mpc consistent with a cos- 
mology with Hq = 50kms _1 Mpc _1 . 
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Residual- Velocity-Regions (HRVR; residual veloc- 
ity field and the shell/chimney candidates, regions 
which have a systematic deviation from rotational 
velocity) were found also in the ionized gas and in- 
terpreted as regions with vertical motions related 
to ongoing star forming processes in the disc on the 
basis of Fabry-Perot Ha observations (Jimenez- 
Vicente & Battancr 2000). Schulman ct al. (1996) 
also computed Lfir = (5.8 ± 1.2) x 1O 9 L and 
L Ha = 1.0 ±0.3 x 10 8 L Q . 

Such high FIR and Ha luminosities (and SFRs 
therefore) would naturally result in a high Super- 
nova rate in this galaxy. Indeed, NGC 5668 is 
known as a "SNe factory" due to the discovery 
of multiple supernova explosions in recent epoch, 
namely SN2004G (Nakano et al. 2004), SN1952G, 
SN1954B (Boffi ct al. 1999). NGC 5668 has also 
been recently observed by a number of instruments 
and facilities, including SAURON at the WHT, 
SDSS, Spitzer and as a part of the Medium-deep 
Imaging Survey of GALEX. This data-set in com- 
bination with the PPAK mosaic obtained as part 
of this work should allow a very detailed analysis 
of the evolution- and chemical-enrichment history 
of NGC 5668. 

3. Observations and Data Reduction 
3.1. IFS Observations 

We have observed the nearly face-on spiral 
galaxy NGC 5668 with the PPAK IFU of PMAS 
at the Calar Alto (CAHA, Spain) observatory 3.5 
m telescope (Kclz et al. 2006). The observations 
were carried out on June 22-24, 2007. We used the 
PPAK mode that yields a total field-of view (FoV) 
of 74" x 65" (hexagonal packed) for each pointing. 
We covered a total area of roughly 2x3 arcmin 2 
with a mosaic of 6 PPAK pointings (see Fig. 1). 
We used the V300 grating covering a wavelength 
range of 3700 -7100 A, (lOAFWHM, correspond- 
ing to a ^300 kms -1 at H/3). With this spec- 
tral configuration we covered all the optical strong 
emission lines in which we are interested. A total 
of 112 individual images were taken during the ob- 
serving run. These images include 19 frames on- 
source for the six pointings obtained in NGC 5668 
(see Table 2 for a summary of the science observa- 
tions), bias frames (22 in total), sky-flats (19), fo- 
cus images (6), HgCdHe arcs (14), tungsten-lamp 
flats (14), observations of the spectro-photometric 



standard star Hz44 (3), dome-flats (6). The spa- 
tial sampling is determined by a hexagonal array 
of 331 densely packed optical fibers for the object 
(science- fibers) , 36 fibers for the sky and 15 cali- 
bration fibers. Each fiber has a diameter of 2.68" 
and a pitch 3.458" with a high filling factor in one 
single pointing (65%). For a more exhaustive de- 
scription of the instrument see Kclz et al. (2006). 

3.2. IFS data reduction 

The reduction procedure applied to NGC 5668 
follows the techniques described in Sanchez (2006) 
using R3D. This is a software package developed 
specifically for the reduction of fiber-based IFS 
data, which has been extensively used for the 
reduction of PMAS data, in combination with 
IRAF 3 packages and Euro3D software (Sanchez 
2004). Six pointings were observed for NCG 5668 
leading to total covered area of 2 x 3 arcmin 2 . 
Each pointing results from the combination of 
three/four images, each one with a exposure time 
of 1000 seconds. Our combined raw mosaic in- 
cludes a total of 2292 spectra (1982+310 sci- 
ence+internal calibration spectra) covering the 
spectral range 3700- 7131 A. The pre-reduction of 
the IFS data consists of all standard corrections 
applied to the CCD that are common to the re- 
duction of any CCD-based data: 

• Bias subtraction: A master bias frame was 
created by averaging all the bias frames ob- 
served during the night and subtracted from 
the science frames. 

• Cosmic-ray rejection: We median combined 
different exposures (at least 3) of the same 
pointing. This procedure is able to remove 
95% of all cosmic-ray events from each point- 
ing. 

After these corrections, the IFS data reduction 
consists of: 

• Spectra extraction: We identified the posi- 
tion of the spectra on the detector for each 
pixel along the dispersion axis. We then co- 
added the flux within an aperture of a 5 pix- 
els around the trace of the spectra and we 

3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 
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extracted the flux corresponding to the dif- 
ferent spectra at each pixel along the disper- 
sion axis. Finally the ID extracted spectra 
are stored in a row-staked-spectra file (RSS). 

• Wavelength calibration: We corrected for all 
distortions and determined the wavelength 
solution by using an arc calibration lamp ex- 
posure. Calibration data cubes were taken 
at the beginning and at the end of the 
night by illuminating the instrument with a 
HgCdHe lamp. The V300 wavelength cal- 
ibration was performed using 14 lines in 
the considered spectral range. The rms 
of the best-fitting polynomial (order 4) was 
0.45 A (i.e. FWHM/20). 

• Fiber-to-fiber transmission: We corrected 
the differences in fibcr-to-fibcr transmission 
by using twilight sky exposures from which 
a Fiber-Flat frame was computed. 

• Sky emission subtraction: For each pointing, 
a sky image of 300 seconds (shifted +5' in 
Dec) was taken. We median combined the 
different sky images to create a sky frame 
which was then subtracted from the science 
data frame. 

• Relative flux calibration: We relative-flux 
calibrated our data using observations of 
the spectrophotometric standard star Hz44 
(a(2000)= 13h 23m 35.37s, <5(2000)= +36d 
08' 00.0", Oke 1990) obtained during the 
night with the same instrumental setup as 
for the object. We applied all the previous 
steps to the calibration star frame and ex- 
tracted the standard star spectrum from the 
brightest 2. 7"- wide PPAK fiber centered on 
the star. Finally, we derived the ratio be- 
tween counts per second and flux. In this 
way we obtained the instrumental response 
which is applied to the science frames in or- 
der to flux-calibrate them. Note that al- 
though this procedure also provides an abso- 
lute calibration for the data this is uncertain 
due to potential flux loses during the obser- 
vation of the standard star and variations in 
the transparency throughout the night. 

• Differential Atmospheric Refraction (DAR) 
correction: We calculated a total theoretical 



DAR value using a Calar Alto public tool 
obtaining a value of 1" (the mean value of 
airmass used is 1.3). In this case the differen- 
tial atmospheric refraction can be considered 
negligible compared with the fiber size (2.7" 
in diameter). Therefore, we considered the 
DAR correction is not required for this data 
cube. In any case we never obtained any flux 
measurements from regions with less than 5 
adjacent fibers in size. 

• Locating the spectra in the sky: A position 
table that relates each spectrum to a certain 
fiber gives the location of the spectra in the 
sky. 

• Mosaic reconstruction: We built a single 
RSS file for the whole mosaic by adding 
pointings previously normalized to a ref- 
erence ('master') pointing. The 'master' 
pointing is that having the best sky sub- 
traction and most optimal observing condi- 
tions regardless of the geometric position of 
the pointing in the mosaic. In our case this 
'master' pointing was number 5 (offsets (0,- 
60), see Fig. 1). This procedure makes use 
of the overlapping fibers between different 
pointings (11 fibers in total) to re-scale each 
new added pointing using the flux ratio in 
the continuum region 5869-5893 A between 
the two pointings (the newly added and the 
'master' one). 

• Absolute Flux Re- Calibration: We also com- 
puted an absolute flux re-calibration based 
on SDSS photometry. Only two of the five 
filters of SDSS ( g and r, X eff = 4694 and 
6178 A, respectively) are used because these 
are the ones for which their pass-bands are 
fully covered by our spectra. We convolved 
the whole datacube of NGC 5668 with the 
SDSS g- and r-filtcr pass-bands to derive 
the absolute spectrophotometry. Firstly, we 
measured the AB magnitude within the foot- 
print of the entire PPAK mosaic and in 7 
concectric annuli for both PPAK and SDSS 
images. Then, these magnitudes were con- 
verted to fluxes using the prescription in the 
SDSS documentation 4 in order to calculate 
the flux ratio in both bands. The resulting 

4 http: //www. sdss3.org/dr8/algorithms/fluxcal.php 
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average scaling factor found is 0.9 for the 
two data pairs. This result implies a factor 
0.04 in magnitudes. Note that the error in 
SDSS photometry alone is estimated to be 
2% mag. In our case we have an uncertain 
of 5% in fluxes so the total absolute error 
associated to line fluxes is 8%. Obviously, 
relative measurements are more precise. All 
results showed in the following analysis (and 
in Table 3) are recalibrated based on these 
SDSS photometry measurements. 

3.3. Imaging observations 

In order to complement the IFS data of 
NGC 5668 presented here, we have also com- 
piled multi-wavelength broadband images ac- 
quired with different facilities: UV imaging data 
from GALEX, optical data from SDSS and near- 
IR data from Spitzer (see Fig. 2). The GALEX 
telescope (Martin et al. 2005) observed NGC 5668 
as part of its Medium-deep Imaging Survey. Im- 
ages at both the FUV (A e/ / = 151.6 nm) and NUV 
(A e // = 226.7 nm) bands are available for this 
galaxy. The final pixel scale is 1.5"/pixel, and the 
PSF has a FWHM of- 5-6" in both bands (Mor- 
rissey et al. 2007). At the distance of NGC 5668, 
this corresponds to a physical scale of —700pc. 
The uncertainty of the zero-point, which is cal- 
ibrated against white dwarfs, is estimated to be 
less than 0.15 mag (Gil de Paz et al. 2007). 

Optical images in the ugriz bands were re- 
trieved from the SDSS (York et al. 2000) Data 
Release 7 archive (Abazajian et al. 2009). The two 
different scans on which NGC 5668 lies were flux- 
matched and mosaicked together using a custom- 
built task in IRAF. We applied the photometric 
calibration specified in the SDSS DR7 Flux Cali- 
bration Guide 5 , relying on the calibration factors 
of the scan used as reference when creating the 
mosaic. According to the information provided 
by the SDSS project, the photometric relative er- 
rors are of the order of 2%-3%. The plate-scale 
is 0.396" per pixel with a PSF FWHM of 1.4" 
(170 pc at the distance of the galaxy). 

Near-IR images at 3.6 ^m and 4.5 /im were 
taken by Spitzer (Werner et al. 2004) using the 
IRAC instrument (Fazio et al. 2004). The cor- 
responding Post Basic Calibrated Data (PBCD) 



5 http://www.sdss.org/dr7/algorithms/fluxcal.html 



were downloaded from the Spitzer archive . These 
images are already flux-calibrated and delivered 
in units of MJysr -1 . The pixel scale is 0.6" per 
pixel, and the PSF has a FWHM of 1.7" in both 
channels, probing a spatial extent of 200 pc (Reach 
et al. 2005). The photometric error (—2%) is dom- 
inated by the aperture corrections that need to be 
applied to account for the scattered light on the 
array. This total uncertainty is estimated to be 
around 10% 7 in this case. 

4. Analysis 

4.1. Spectroscopy 

Many important topics in astrophysics involve 
the physics of ionized gases and the interpretation 
of their emission line spectra. Powerful constraints 
on theories of galactic chemical evolution and on 
the star formation histories of galaxies can be de- 
rived from the accurate determination of chemi- 
cal abundances either in individual star forming 
regions or distributed across galaxies or even in 
galaxies as a whole. In this sense chemical abun- 
dance are a fossil record of its star formation his- 
tory. In addition, the distribution of H II regions 
is an excellent tracer of recent massive star forma- 
tion in spiral galaxies. The previous evolutionary 
history of the gas, and the local ionization, den- 
sity and temperature are important features to un- 
derstand the physical conditions prevailing in the 
emitting regions where they are emitted. 

For all these reasons, our first aim is to iden- 
tify H II regions in our target, or at least H II 
complexes (see Fig. 3) 8 . We extract from the 
PPAK data cube the emission line fluxes corre- 
sponding to each of these regions. Furthermore, 
we also compute 2D maps of the flux of strong 
emission lines, such as [O n] AA 3726,3729AA 9 ; 
[O in] AA 4959,500711; H/3; [Nil] A 65481; Ha; 
[Nil] A 65831 and [Snj AA6717,673lll. 



6 Spitzer Proposal ID 69: A Mid-IR Hubble Atlas of Galax- 
ies, PI: G. Fazio 

7 http://ssc.spitzer.caltech.edu/irac/calib/extcal/ 
8 Hereafter we use the terms 'H II regions' and 'H II com- 
plexes' indistinctly, even although individual H II re- 
gions are not spatially resolved at the physical resolution 
achieved in our IFS data. 

9 Thc [On] AA 3726,3729AA doublet is spectrally unresolved 
in our data, therefore and hereafter we will refer to it as 
[O n] A3727A. 
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Fig. 2. — Left: IRAC color image of NGC5668, using the 3.6/im and 4.5^tm images. Since these bands primarily probe old 
stars, no significant color variations are found. The red dashed lines show the width of the weak bar visible in the optical images 
(~12" in size) and the blues dashed circle corresponds to the break radius of the metallicity gradient (r~36"). Right: GALEX 
color image, using the FUV and NUV bands as the blue and red channels, respectively, and a linear combination of both in the 
green one. The diameter of the yellow dashed circle is the optical one quoted in the RC3 catalog, at 25.5 mag/arcsec 2 in the B 
band (~3.3'). Both frames are 250" in size. North is up and East is to the left. 



Jf.. 1.1. Region selection 

In spiral galaxies, H II regions are strongly con- 
centrated along the spiral arms and in the galac- 
tic plane and are the best objects for tracing the 
structure of spiral arms in galaxies. We are in- 
terested in the study of the optical emission lines 
from H II regions because they represent the pri- 
mary mean of performing gas-phase diagnostics in 
galaxies (Pagel 1986, Osterbrock & Ferland 2006). 
These data were analyzed using IRAF and E3D 
(Sanchez 2004) tasks and some custom-made IDL 
(Interactive Data Language) scripts. In order to 
select the H II regions of NGC 5668 we com- 
puted a synthetic Ha image (the brightest emis- 
sion line in the optical spectrum of H II regions un- 
der most physical conditions) from the data cube 
using E3D. The spectral window for the Ha line 
was defined as a narrow wavelength range from 
6589 A to 6613 A. Two bands (6457-6553 A and 
6667-6730 A) close to Ha line were also extracted 
for continuum subtraction. We normalized the 
"continuum" maps to match the width of the Ha 
intensity map. From the synthetic Ha image we 



can then identify H II regions. We visually se- 
lected 62 H II complexes in the galaxy from this 
image (see Fig. 3). The apertures were defined so 
they would approximately reach the same surface 
brightness level in Ha. We also imposed that these 
apertures would be large enough to include at least 
5 PPAK fibers (in order to minimize atmospheric- 
refraction effects) and that each region would be 
isolated from bright neighbors. As a result of the 
conditions imposed the majority of regions are the 
sum of 5 PPAK fibers, except in the case of regions 
number 17, 37 and 62, which are the superposition 
of one bright and several, marginally-resolved faint 
H II regions, and region 58 which extends 15" in 
the SE-NW direction approximately following the 
spiral pattern of the galaxy northern arm. See 
Fig. 3 for the positions and sizes of the regions of 
interest selected for this study. We extracted the 
spectrum coming from each region by averaging 
the signal for each fiber belonging to a particular 
region. 

In many spiral galaxies of early and intermedi- 
ate Hubble type (Sa-Sc), active star formation is 
organized in a ring-like structure that often con- 
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tains a large fraction of the entire star formation 
activity of the galaxy. To investigate this and the 
ubiquitous radial variation of the physical prop- 
erties in spiral galaxies, we also selected 18 con- 
centric annuli, centered on the peak of the optical 
continuum emission. A width of about 5" proved 
to be a good compromise in term of the trade-off 
between spatial resolution and depth. The outer- 
most ring analyzed is located at i?i as t = 95". 

In a first step we analyze the data to verify the 
quality of the spectra extracted in all regions: we 
identify and eliminate bad fibers (e.g. those placed 
on bright field stars or affected by cross-talk), we 
interpolate over the spectrum at the position of 
the brightest sky lines to clean out our spectra. 
We trim the wavelength range in each spectrum to 
avoid the poor instrumental efficiency at the very 
edges. This leads to a useful wavelength range of 
3702-6997 A. 

4-1.2. Emission-line fluxes 

The radiation emitted by each element of 
volume in a region depends upon the abun- 
dance of the elements, determined by the pre- 
vious evolutionary history of the galaxy, and 
the local ionization, density and temperature. 
The most prominent spectral features are emis- 
sion lines, many of these are collisionally-excited 
lines. In the V300 grating setup the main emis- 
sion lines observed are: [On]A3727A doublet; 
[Om] AA4959,5007AA; H/3; [Nn]A6548A; Ha; 
[Nn]A6583A and [S n] AA 6717,6731 AA. 

We identified two different types of regions ac- 
cording to the intensity of the continuum: (1) 
regions where the continuum was bright enough 
so some absorption features were clearly visible 
and, therefore, the spectrum of the underlying 
stellar population could be fit using stellar popula- 
tions synthesis models and (2) regions with fainter 
and/or noisier continua where no spectral absorp- 
tion features could be identified. 10 

We refer to section 4.1.3 for a description of the 
procedure followed to fit the underlying continuum 
in the case of the former regions; here we discuss 



The limit in the continuum intensity between the two types 
of spectra was set to 1 X 10 -16 ergsec -1 A~ 1 . For those 
spectra where the continuum was fainter (in either H/3 
or Ha) than this number no attempt to carry out a full- 
spectral fitting to the continuum was attempted. 



the procedure followed to derive accurate emission 
line fluxes in regions with negligible or noisy con- 
tinua. Initially, we have identified 73 H II regions 
in our synthetic Ha image but we excluded some 
of these regions because of their low SNR or be- 
cause they were later identified as bad columns or 
field stars. 

We first masked out the main emission lines 
and used a third-order polynomial function to 
fit the continuum. We then selected two spec- 
tral ranges excluding emission lines in three dif- 
ferent wavelength ranges [On] A3727A; H/3; Ha 
regions and measured the continuum level in each 
of these ranges. In a second step we averaged 
the two continua and subtracted the result from 
each emission line spectra in order to obtain a 
decontaminated flux. We analyzed the spectra 
using the IRAF package ONEDSPEC and STSDAS 
packages. For well isolated emission lines, such 
as [On]A3727A, US, H 7 n , [Om] A5007A, we fit 
the emission line profile to a single gaussian func- 
tion. For partially-blended emission lines, such as 
the [Nn]A6548A; Ha; [Nn]A6583A triplet, and 
[S n] AA 6717,673lAA doublet, we used the task 
ngaussfit that fits simultaneously multiple gaus- 
sian functions. We calculated the typical signal- 
to-noise ratio (SNR) of our data. In the case of 
[O n] A 3727A the SNR ranges between 3.4 and 23. 
In the case of H/3 we obtain similar SNR values, 
ranging between 2.5 and 23. For the part of our 
spectra where Ha emission line is located the SNR 
values varies between 6 to 67. 

Before analyzing the data in terms of chemi- 
cal abundances or star formation rates, the ob- 
served emission-line fluxes must be corrected for 
various effects. For example, the presence of dust 
between the zone of emission and the observer or 
the possible underlying stellar absorption in the 
case of the hydrogen recombination lines can al- 
ter considerably the intrinsic emission line fluxes. 
We correct for underlying stellar absorption in the 
hydrogen Balmer lines using the values obtained 
from the analysis of the ring spectra (see section 
4.1.3 for fitting details). These best- fitting equiv- 
alent widths in absorption show very little varia- 
tion with radius ranging between -1.9 and -2.3 A 
in H/3 and -1.7 and -1.8 A in Ha. We apply an av- 



Note that we detect these lines only in 23 HII Regions, so 
we not use H<5 and H7 to derive physical properties. 
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Fig. 3. — PPAK continuum-subtracted Ha emission map. The H II regions selected in NGC 5668 are shown as blue circles 
see section 4.1.1 for details. The red dashed lines indicate the radii of the rings and the red cross represents the central fiber of 
the mosaic. North is up and East is to the left. 
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erage correction of -2 and -1.7 A, respectively for 
H/3 and Ha, to the spectra of the 58 H II regions 
where the continuum emission was too faint for 
carrying out a full spectral fitting. The emission 
fluxes were then corrected for reddening using the 
Ha/H/3 Balmer decrements after adopting an in- 
trinsic ratio of 2.86 (Osterbrock & Ferland 2006). 
Additionally, we verified that in the case of regions 
of high equivalent width in emission (typically no 
less than 5 A in H/3) the corrected H/3/H7 ratios 
were consistent with the predictions for case-B re- 
combination value at a typical T e of ^10 4 K. In 
Table 3 we present the results obtained from the 
analysis of H II regions and concentric annuli in 
NGC 5668. 

Following a similar analysis on a pixel-by-pixcl 
basis, we also generate maps of ionized-gas extinc- 
tion, radial velocity, emission line fluxes and stellar 
absorption equivalent width. In Fig. 4 we present 
the maps of the Ha line flux; [0 11] A 3727A dou- 
blet flux; Ha continuum intensity and equivalent 
width (EW). Note that while the EW(Ha) in emis- 
sion gets locally higher as we move towards the 
outer parts of the disk, this is mainly due to the 
decrease in the intensity of the adjacent underly- 
ing (i?-band) continuum. Indeed, when we com- 
pute the azymuthally-averaged EW(Ha) from the 
ring spectra this changes less dramatically, with 
values ranging from 20 to 50 A in all cases ex- 
cept for the very outer rings (see Table 4). The 
very high EW(Ha) values measured in the rings 
beyond ~ 70 " are likely due to strong Ha emis- 
sion associated to the bright H II complex located 
in the galaxy Northern spiral arm. Note that in 
the very outer disk, where the total SFR is low, 
stochasticity in the number and luminosity of H II 
regions might lead to significant fluctuations in 
the EW(Ha) compared with values averaged over 
timescales of a few hundred Myr or even a rotation 
period. 

4-1.3. Stellar-absorption correction 

As discussed above, for the rings and some H II 
regions the continuum is sufficiently strong that 
allows us to clearly identify absorption features 
in their spectra. In this way we can fit the ab- 
sorption and emission components in the spectra 
simultaneously. In this case we made use of a own 
IDL script that finds the best-fitting stellar popu- 
lation synthesis model among the ones in libraries 



by Bruzual & Chariot (2003) (hereafter BC03) 
and Sanchez-Blazquez et al. (2006) (also known 
as MILES). This procedure allowed us to mea- 
sure intensities for the emission lines already cor- 
rected for stellar continuum absorption (nebular 
continuum emission is negligible at the EW(Ha) 
measured anywhere in this object). Both libraries 
yield good fits to our spectra, yet the analysis pre- 
sented here is based on the results obtained us- 
ing the BC03 library since this provides a good 
coverage of physical parameters and good spectral 
resolution and has been more extensively tested 
than MILES. We explore the whole range in age 
and metallicity in the BC03 library for simple stel- 
lar populations (SSPs) with a spectral resolution 
of 3 A across the whole wavelength range from 
3200 A to 9500 A. These templates have 11 val- 
ues in age (5, 25, 100, 290, 640, 900 Myr, and 1.4, 
2.5, 11, 13, 17 Gyr) and for each of them six dif- 
ferent metallicity values: z=0.0001, 0.0004, 0.004, 
0.008, 0.02, 0.05. Our fitting program is able to 
find the model that best fits the underlying stellar 
population (stellar continuum) in this BC03 stel- 
lar library by minimizing the residuals between 
the model and observed spectra for a set of input 
parameters. The program finds the best fitting pa- 
rameters for the recession velocity, velocity disper- 
sion, normalization factor and stellar-continuum 
attenuation in the V band for each region. 

In order to attenuate the stellar continuum as 
a function of wavelength we adopt a dependence 
of the optical depth with wavelength of the form 
A~ - 7 , which, according to the results of Chariot 
& Fall (2000), is a good approximation to the ex- 
tinction law in star forming galaxies at these wave- 
lengths. This analysis provides us with the best- 
fitting SSP model for each spectrum along with 
best-fitting values for the recession (tropocentric) 
velocity, velocity dispersion, normalization factor, 
and continuum dust attenuation. 

Note that the age and metallicities derived are 
in general highly uncertain but, on the other hand, 
the Ha and H/3 equivalent widths in absorption 
are much more precise (see Marmol-Qucralto et al. 
2011 for a detailed description on the feasibility 
in the derivation of the properties of the under- 
lying stellar population). The residual spectrum 
obtained after the subtraction of the best-fitting 
SSP model is then used for the subsequent analy- 
sis of the emission-line fluxes. 
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Fig. 4. — NGC 5668 synthetic maps. Top: Emission lines maps at Ha (left) and [On] A3727A doublet (right). Bottom: 
Ha continuum (left) and EW (right); we calculated these maps assuming an average absorption equivalent width in H« of -2 A. 
Ha emission line fluxes in this plot are represented by the isocontours. The mosaic center is marked with a cross. North is up 
and East is to the left in all cases. 
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Fig. 5. — PPAK spectra extracted of the rings selected in NGC 5668. From bottom to top we show the spectra of rings 
with increasing radius. The spectra are obtained summing all fibers within each ring. Two bright sky-lines residual ranges are 
masked. These spectra are plotted in linear scale and shifted in the y-axis for a better visualization. 
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We calculate the errors in all line fluxes from 
the rms in the spectra (in the region adjacent to 
each emission feature) after the best-fitting stellar 
synthesis model has been subtracted. We assume 
that these rms measurements were due mainly to 
photon noise and to sky subtraction. These errors 
were scaled to those spectra for which the fainter 
continuum emission prevented earring out a full 
spectral fitting. Finally, we used the errors in the 
continuum to calculate the errors associated to the 
line fluxes 12 . 

4.2. Integrated spectrum 

Some advantages of using IFUs are that all 
the spectra are obtained simultaneously and that 
we can use the IFU as a large-aperture spectro- 
graph. From our datacubc we add up the spec- 
tra of all the fibers to create an integrated spec- 
trum of NGC 5668 (see Fig. 6). In this figure we 
also compare this spectrum with the correspond- 
ing fluxes derived from the Sloan Digital Sky Sur- 
vey photometry. Note that the g' and r' bands 
overlap completely with our spectral range. Since 
the spectroscopy data were not taken under photo- 
metric conditions and given that some flux losses 
are expected as part of the PPAK observations of 
our spcctrophomctric standard star, a 8% differ- 
ence was found between the two datasets. This 
offset, which is wavelength independent, has been 
already applied to the spectrum shown in Fig. 6 al- 
though was not applied to our spectral datacube 
and the line fluxes measured hereafter. Similar 
wavelength-independent offsets are also found in 
other works based on PPAK spectroscopic obser- 
vations (Rosales-Ortega 2009). In any case, this 



The errors are calculated from the residual fit-models of the 
underlying stellar population; we are assuming that these 
residuals are primarily due to photon noise and sky subtrac- 
tion. We note that we have pointings with a slightly differ- 
ent spatial resolution and several combinations of exposure 
time and number of individual exposures; for this reason 
the continuum level emission and the associated noise was 
derived for each pointing separately. The sky subtraction 
to the spectra of NGC 5668 was done, see the data reduc- 
tion steps in Sanchez (2006), by median average the 331 
sky-spectra. The resulting spectrum was then subtracted 
from the cube of NGC 5668. The analysis of these sky- 
subtracted spectra showed that the error in the continuum 
level emission increases linearly with the number of fibers 
(but not to the square root). This indicates that for aper- 
tures including a large number of fibers the total error is 
dominated by the sky-subtraction error. 



does not have any impact on the physical proper- 
ties derived (dust attenuation, oxygen abundance, 
electron density) as these are determined by the 
ratio of emission-line fluxes alone. In Table 3 (last 
row), we present the emission-line flux ratios found 
for the integrated fibers spectrum of NGC 5668. 

4.3. Imaging 

Prior to the measurement of the surface bright- 
ness profiles in NGC 5668, foreground stars, back- 
ground galaxies and cosmetic artifacts had to be 
masked out from our multi-wavelength imaging 
data. Radial profiles were obtained using the 
IRAF task ellipse. We measured the mean in- 
tensity within concentric circular annuli centered 
on the galaxy's nucleus, using a radial increment 
of 6" between adjacent annuli. Photometric errors 
were derived as explained in Gil de Paz & Madorc 
(2005) and Munoz-Mateos et al. (2009a). The un- 
certainty of the mean intensity at a given radius is 
computed as the quadratic sum of two terms: the 
poissonian noise of the galaxy's light and the er- 
ror in the sky level. The latter being the result of 
the combination of both pixel-to-pixel variations 
an large-scale background errors. 

5. Results 

5.1. Attenuation 

More than 80 years have passed since Trum- 
pler's discovery of color excesses provided the first 
definitive proof of the existence of interstellar dust 
(Trumpler 1930). However the nature of interstel- 
lar dust still remains unclear and how dust both 
reddens and attenuates the light from the stars is 
one of the least understood of the physical phe- 
nomena which take place in galaxies. 

The dust content of galaxies is a critical is- 
sue given its important impact on the observa- 
tional properties of galaxies (mainly in the opti- 
cal and UV). The extinction in a galaxy depends 
first on the amount of dust and its composition, 
and also on the distribution of dust relative to the 
light sources (Calzetti 2001). For instance, some 
authors favor a foreground screen dust geometry 
model (Calzetti ct al. 1994, Calzetti et al. 1996), 
while others propose hybrid models with the dust 
partially distributed in a foreground screen and 
partially concentrated in the star-forming regions 
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Fig. 6. — The integrated spectrum of NGC 5668 is shown as a black line. The SDSS g' and r' band photometry data are 
shown as grey squares. Horizontal error bars represent the FWHM of each filter. 
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(Chariot & Fall 2000). In this subsection we dis- 
cuss the methods used to estimate dust attenu- 
ation. We calculate the attenuation of the stel- 
lar continuum from the UV data available on 
NGC 5668 and these estimates are then compared 
with the ionized-gas attenuation derived from the 
Balmer decrement in both individual H II regions 
and concentric annuli. 

5.1.1. UV- continuum attenuation 

In this case we estimate the attenuation of 
the stellar continuum using UV data alone. Be- 
cause the UV radiation is preferentially emitted 
by young stars (~100Myr) and because the dust is 
most efficient in attenuating UV light, rest-frame 
UV observations can lead to incomplete and/or 
biased reconstructions of the recent star forma- 
tion activity and star formation history of galax- 
ies where dust absorption is expected to be signif- 
icant, unless proper corrections are applied. Ra- 
diative transfer models suggest that the total-IR 
(TIR, 3-1100/xm) to UV luminosity ratio method 
(i.e. Buat 1992; Xu & Buat 1995; Meurer et al. 
1995; Meurer et al. 1999) is the most reliable es- 
timator of the dust attenuation in star-forming 
galaxies because it is almost completely indepen- 
dent of the extinction mechanisms (i.e. dust/star 
geometry, extinction law, see Xu & Buat 1995; 
Meurer et al. 1999; Gordon et al. 2000; Witt & 
Gordon 2000). However, this would require hav- 
ing both far-infrared data at the same resolution 
as the UV data. In the case of galaxies at the dis- 
tance of NGC 5668 this would provide very little 
spatial information on the radial variation of the 
dust attenuation even if data from state-of-the- 
art infrared facilities such as Herschel would be 
available. Even though IRAC and MIPS images 
for NGC 5668 are available in the Spitzer archive, 
so theoretically, we can obtain profiles in 8, 24, 
70 and 160 ^m but the poor resolution at 160 /jm 
(38") would yield a total-IR profile with just a cou- 
ple of data-points. Under certain circumstances, 
the UV attenuation can be indirectly estimated us- 
ing the slope of the UV spectrum (denoted as /3), 
in the sense that redder UV colors are indicative of 
larger attenuations. This so-called IRX-/3 relation 
has been widely used in starburst galaxies, where 
most of the UV light comes from newly-born stars 
(Calzetti et al. 1994; Heckman et al. 1995; Meurer 
et al. 1999). In normal star- forming spirals, more 



evolved stars can also contribute to the UV flux, 
shifting the IRX-/3 relation to redder UV colors 
and increasing the overall scatter (Bell 2002; Buat 
et al. 2005; Seibert et al. 2005; Cortcsc et al. 2006; 
Gil de Paz et al. 2007; Dale et al. 2007). The use of 
radial profiles instead of integrated measurements 
seems to reduce the global scatter (Boissier et al. 
2007; Munoz-Mateos et al. 2009b). 

Here we apply the IRX-/3 relation provided by 
Munoz-Mateos et al. (2009b) to FUV-NUV mea- 
surements obatined from the GALEX images of 
NGC 5668. This relation was specifically cali- 
brated for normal nearby spirals, so it should be 
applicable to NGC 5668 as well. However, before 
applying this technique to subregions, we checked 
this assumption using the total galaxy IR lumi- 
nosity of NGC 5668 from the IRAS fluxes quoted 
by Wang & Rowan-Robinson (2009) leading to 
a total IR flux of log(F T i R ) = 13.64 x 10~ 23 
(ergs -1 cm -2 ). This value combined with the 
FUV flux obtained from our GALEX FUV im- 
age leads to log(TIR/FUV) - 0.18. On the other 
hand, for a global UV color of FUV-NUV=0.27, 
the IRX-/3 yields log(TIR/FUV) - 0.39. This is 
somewhat larger than the observed value, but still 
lies within the la scatter around the IRX-/3 rela- 
tion of Munoz-Mateos et al. (2009a). Thus, us- 
ing the radial profiles described in Section 4.3, we 
determined the radial variation of the FUV-NUV 
color in steps of 6". The UV color was then used 
to obtain profiles of TIR/FUV and TIR/NUV ra- 
tios using the prescriptions of Munoz-Mateos et al. 
(2009a). These values were then transformed into 
dust attenuation in the FUV (Apuv) and NUV 
(Ajvt/v) by means of the recipes by Cortese et al. 
(2008), which take into account the IR excess as- 
sociated to the extra dust heating due to the light 
from old stars. The resulting attenuation profiles 
are shown in Table 5. Note that the absolute un- 
certainty in the TIR/FUV ratio from the FUV- 
NUV is larger than the relative change in such 
ratio across the galaxy, see Munoz-Mateos et al. 
(2009a). 

Now that we have derived the extinction in 
the UV, we can calculate the extinction in other 
bands using an attenuation curve. Such redden- 
ing corrections have normally been undertaken 
using a number of extinction curves, including 
those of Seaton (1979) (in the UV), Cardelli et al. 
(1989), Savage & Mathis (1979), Ardeberg k Vird- 
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efors (1982) and Fitzpatrick (1999). These "stan- 
dard" curves are often used interchangeably, on 
the understanding that they should give broadly 
similar results. The parameterized curve by 
Cardelli, Clayton, & Mathis (1989) law is com- 
monly used to fit the extinction data both for dif- 
fuse and dense interstellar medium. In our case, 
the continuum attenuation is obtained using at dif- 
ferent wavelengths two different extinction laws, 
Cardelli et al. (1989) and Calzetti (2001). We 
based our calculation on the NUV data, adopting 
the relation with the color excess E(B-V) given by 
Cardelli et al. (1989), E(B-V)= A NUV /8.0 and for 
Calzetti et al. (1994), E(B-V)= A OT y/8.22 (Gil 
de Paz et al. 2007). We then adopt the average 
values of Rv.cardeiu =3.1 and Rv,caizetu =4.05. 
For sake of clarity Fig. 7 plots only the results ob- 
tained in the case of the Calzetti extinction law 
as the two curves are consistent within the errors 
(but offset by -0.09 magnitudes in the case of the 
Cardelli law). 

5.1.2. Ionized-gas emission-lines attenuation 

We now made use of the Balmer decrement in 
individual H II regions and concentric annuli for 
determining the amount of dust extinction in the 
ionized-gas emission along the line of sight. The 
relative ratios of the Balmer lines of hydrogen are 
often used as extinction indicators due to the fact 
that they are observationally convenient (being in 
the optical band), strong and their intrinsic rel- 
ative flux ratios are fairly well determined from 
atomic theory. Under Case B conditions with tem- 
peratures ^10,000 K and electron densities ~100 
cm 3 , the theoretical Ha/H/3 ratio should be close 
to 2.86 (Osterbrock & Ferland 2006). Since extinc- 
tion is more severe in H/3 than at Ha wavelengths, 
the net effect is to increase the observed Balmer- 
line ratio. As we discussed above, the results are 
dependent of the dust model of choise. Here we as- 
sume that the dust is located in a uniform screen 
between us and the gas and that the extinction 
law in the optical is similar to the Galactic one, as 
parametrized by Cardelli et al. (1989). 

5.1.3. Optical- continuum attenuation 

The last method that we have used for the 
extinction is based in our IDL fitting-program 
for the underlying stellar population (see 4.1). 



From the output values for each concentric an- 
nuli we find a mean value of Ay^l.7 mag. This 
value has a very large uncertainty due to the lim- 
ited wavelength range of our spectroscopic ob- 
servations which leads to strong age-metallicity- 
extinction-SFH degeneracies (see Gil de Paz & 
Madore 2002). The results are summarized in 
Fig. 7. In general, the ionized-gas dust attenu- 
ation has a mean value of ~1.5 magnitude, which 
is in agreement to what is found in other spi- 
ral disks (Gil de Paz et al. 2007). We find that 
the gas attenuation is larger than the continuum 
attenuation by <E(B-V) con ti„ Mttm > =0.36 x < 
E(B-V)i on i zec ;- g as > when the Calzetti attenua- 
tion law is adopted; note that ionized-gas color 
excess (derived from the Balmer decrement) is 
corrected for a Milky Way foreground reddening, 
E(B - V) = 0.037. This is somewhat expected 
given that, as shown in Calzetti et al. (2000) and 
Stasihska & Sodre (2001), in nearby spiral galax- 
ies there is evidence that stars, gas and dust are 
typically decoupled (see e.g. Mai'z-Apellaniz et al. 
1998) so the attenuation inferred from the Ha/H/3 
ratio is typically higher than that inferred from 
the spectral continuum of the same wavelength 
(e.g., Calzetti et al. 1994; Mayya & Prabhu 1996; 
Calzetti 1997). Poggianti et al. (1999) and Pog- 
gianti & Wu (2000) explain such results as due 
to selective-dust-extinction effects where a large 
fraction, but not all, of the dust in galaxies is asso- 
ciated with star formation regions, which absorb 
a significant fraction of the light emitted by the 
young stars. 

5.2. Electron density 

The electron density, N e , is one of the key phys- 
ical parameters that characterize gaseous nebulae. 
Most of the density estimates found in the liter- 
ature are based on measurements of a sensitive 
emission-line ratio. In the presence of internal 
variations of electron density, however, these single 
line-ratio measurements may not be representative 
of all ionizing zones. The measurements are based 
on the fact that a relation exists between the colli- 
sional de-excitation of atoms and electron density, 
(Osterbrock & Ferland 2006). If an ion emits sim- 
ilar amounts of energy from different energy levels 
at nearby wavelengths, then the ratio of the emis- 
sion line intensities can be used to obtain an esti- 
mate of the electron density and there are two line 
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Fig. 7. — Left: Representation of the radial variation of dust attenuation. Filled circles represent the values obtained from 
individual H II regions based on the Balmer decrement measurements. The black line represents Balmer decrement data for 
the concentric annuli. The red line is the attenuation (and ±lcr error band) calculated from the slope of the UV continuum (or, 
cquivalcntly, the FUV-NUV color) adopting a Calzetti extinction law. Right: Dust extinction map (MW + internal) of the 
ionized gas in NGC 5668 computed from the Hct/H/3 ratio. The Ha continuum values are represented as isocontours. North is 
up and East is to the left. 



ratios that can be used: [On] AA 3726,3729AA or 
[S 11] AA 6717,673lAA. In our case, we make use of 
the [Sn] AA6717,673lAA line ratio and the elec- 
tron density of the region responsible for the [S 11] 
emission can be determined as: 



I([5ri]A67n) 
fl([5ll]) - J([Sn]A6731) 



1.49 



l + 3.77aA 
1 + 12.8a: J 
(1) 

where x, the density parameter, is defined as 
x = 10~ 4 • n e ■ i" 1 / 2 . Solving this equation for 
the density parameter we can obtain N e by as- 
suming T e — T[5 ii] /10 4 , where Tjsii] is the elec- 
tron temperature of the region responsible for the 
[S11] emission (McCall et al. (1985)). Thus, in 
order to properly estimate this electron density 
a previous knowledge on the electron tempera- 
ture of the region responsible for the [S 11] emis- 
sion lines is required. Assuming the calibration 
of T[NII] as a function of the R23 line ratio (see 
equation 2) given by Thurston et al. (1996) and 
an average difference of 3000 K between T[Nn] 
and T[Sil] (with the former being lower; Garnett 
(1992)) we estimated T[ S ni in the equation above. 
Note that although this estimate for the tempera- 
ture is not very precise (as it does not rely on the 



use of temperature-sensitive line ratios) is accurate 
enough for correcting the densities obtained from 
the [S11] AA6717,673lAA line ratio for tempera- 
ture effects. The large uncertainties associated to 
the determination of T[Sn] prevent us from ex- 
tracting further conclusions from the density mea- 
surement derived. Despite that, we find a mean 
value for N e of 190 cm~ 3 for the H II regions, 
in agreement with the typical value densities in 
H II regions (N e ~10 2 cm -3 , Osterbrock & Fer- 
land 2006). 

5.3. Chemical abundances 

H II regions identify the sites of recent mas- 
sive star formation in galaxies. The rapid evo- 
lution of these stars, ending in supernovae ex- 
plosions, and the subsequent recycling of nucle- 
osynthesis products into the interstellar medium, 
make of H II regions essential probes the present- 
day chemical composition of star-forming galaxies 
across the Universe. The study of nebular abun- 
dances is therefore crucial for understanding the 
chemical evolution of galaxies. Obtaining a di- 
rect measurement of chemical abundances in H II 
regions requires a good estimate of the electron 
temperature. Unfortunately, this implies detect- 
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Fig. 8. — Radial abundance gradient in NGC 5668, where filled-symbols correspond to individual H II regions while the grey- 
shaded area represents the values obtained for the concentric annuli. Top panel: The red line symbolizes the linear regression fit 
to the H II regions data. We obtained a single gradient of 0.035 ± 0.007 dcx/kpc. The blue dashed line represents the gradient 
obtained from the models of Boissier & Prantzos (1999, 2000), with an offset of 1 dex in the metallicities plotted (see text for 
details). Bottom panel: In this plot we draw two types of fits: the grey line reproduces a double fit weighted by the errors; the 
gradient has a value of 0.140 ± 0.016 for inner part and 0.002 ± 0.019 dex for outer part. The grey dashed line reproduces 
the unweighted double fit that yields gradients values of 0.1073 ± 0.0165 (inner) 0.002 ± 0.019 (outer). The color coding of 
the points represents the different line ratios used to calculate the metallicity for the central point of the ratios probability 
distribution (based on Kewley &; Dopita 2002 recipe): green points are calculated from [Nil] /[On], red triangles represent the 
values obtained with [Nil] /Ha and blue diamonds show i?23-based values. Note that the presence of three uncertain points 
(regions 13, 17, 59 have an error major of 0.3) has not influence on the fitting parameters values. 
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ing very faint auroral lines such as [O in] A 4363A 
or [Nn] A5755A, that are very faint at the abun- 
dance levels of most spiral disks. Thus, we must 
use instead calibrations based on the predictions 
of photo-ionization models or on empirical mea- 
surements for strong lines. Among the various 
strong-line methods, the R23 indicator originally 
proposed by Pagcl et al. (1979) stands out as ar- 
guably the most popular. 

_ f([Om]AA4959, 5007) + /([07/]AA3726, 3729) 
23 " JiW) 

(2) 

Like all strong-line diagnostics, R23 as an abun- 
dance indicator has a statistical value, based on 
the fact that the hardness of the ionizing radia- 
tion correlates with metallicity. Numerous cali- 
brations of i?23 in terms of the nebular chemical 
composition can be found in the literature. We 
choose as best set of values for the oxygen abun- 
dance those given by the Kewley & Dopita (2002) 
recipe, which basically adopts a calibration of the 
[N 11] /[0 11] ratio for oxygen abundances above half 
solar and their own calibration of R23 for abun- 
dances below that value 13 . The R23 abundance 
diagnostic depends strongly on the ionization pa- 
rameter, for this reason we also calculated it fol- 
lowing the indications given in Kewley and Do- 
pita (2002). These authors used a combination 
of stellar population synthesis and photoioniza- 
tion models to develop a set of ionization param- 
eter and abundance diagnostics based only on the 
use of strong optical emission lines. These tech- 
niques are applicable to all metallicities. In par- 
ticular, for metallicities above half solar, the ra- 
tio [Nil] /[On] provides a very reliable diagnostic 
since it is ionization-parameter independent and 
does not have a local maximum. This ratio has not 
been used historically because of concerns about 
reddening corrections. However, the use of clas- 
sical reddening curves is sufficient to allow this 
[Nil] /[On] diagnostic to be used with confidence 
as a reliable abundance indicator. Note that the 
calibration of the [Nil] /[On] line ratio as a di- 
agnostic of 12 + log(0/H) relies on the depen- 
dence of the (N/O) abundance ratio on the oxy- 
gen abundance. Thus, no attempt has been made 



We used the value of 12 + log(0/H) Q = 8.93 as in Kewley 
& Dopita (2002). 



in this work to derive the Nitrogen abundance in 
this work. 

The iterative method used as part of this recipe 
allows both these parameters to be obtained with- 
out the need of using temperature-sensitive line 
ratios involving very faint emission lines, that are 
particularly elusive in the case of spiral disks. One 
drawback of using i? 2 3 (and many other emission- 
line abundance diagnostics; see Perez-Montero & 
Diaz 2005) is that it depends also on the ioniza- 
tion parameter (/defined here as q = Suo/n, where 
Shq is the ionizing photon flux through a unit 
area, and n is the local number density of hydrogen 
atoms. Some calibrations have attempted to take 
this into account (eg. McGaugh 1991), but others 
do not (eg. Zaritsky et al. 1994). Another diffi- 
culty in the use of R23 and many other emission- 
line abundance diagnostics is that they are double 
valued in terms of the abundance (12 + log (O / H)). 
Thus, one of the main difficulties in its adoption 
is related to the necessity of locating which of 
two branches (upper and lower) a given H II re- 
gion belongs to, since R23 is degenerate (Bresolin 
et al. 2009a). This is because at low abundance 
the intensity of the forbidden lines scales roughly 
with the chemical abundance while at high abun- 
dance the nebular cooling is dominated by the in- 
frared fine structure lines and the electron tem- 
perature becomes too low to collisionally excite 
the optical forbidden lines. When only double- 
valued diagnostics are available, an iterative ap- 
proach which explicitly solves for the ionization 
parameter, as the one used here that is based on 
the Kewley and Dopita logical flow diagram, helps 
to resolve the abundance ambiguities. In Table 7 
we present the results obtained. The metallicity 
values vary between 12 + log(0/H) = 8.15 and 
12 + log(0/H) = 8.7 (i.e. from approximately 
1/3 the solar oxygen abundance to nearly the so- 
lar value). The lack of some regions in this table 
is due to the poor signal-to-noise or simply not 
detection of H/3 emission (12 regions). 

Simulations of the line flux errors were car- 
ried out in order derive the errors in the oxygen 
abundances. These simulations assumed a gaus- 
sian probability distribution for each line flux and 
that these were not correlated between the differ- 
ent lines. The resulting errors, which are shown in 
Table 7, could then be considered as upper limits 
to the actual line-flux errors (where a contribution 
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of correlated errors is expected). 

In order to quantify the abundance gradi- 
ent, we have carried out three linear regres- 
sions to the data-points, the first one is a fit 
to all points weighted by their errors (xled = 
1.45); the gradient has a value of —0.035 dex/kpc, 
(—0.0042 dex/arcsec). The second and the third 
fits are double fits to the data with a free pa- 
rameter, the radius of break, but one is weighted 
by the errors (black solid line) and the other is 
unweighted (black dashed line; yield both Xred = 
0.54). 

The Kewley & Dopita recipe allow us to make 
use of the most optimal abundance indicator 
among [Nil] /[On], [Nil] /Ha and R 23 . Wc find 
that the three abundance diagnostics (showed 
in the bottom panel of Fig. 8 with different col- 
ors) are homogeneously mixed, Fig. 8. Thus, the 
[N n]/Ha and i?23-based abundance values are well 
mixed and can be derived at almost any galaxy 
radius, while the [Nil] /[On] indicator is lost in 
some H II regions, mostly at large radii; this is, as 
pointed out in Kewley & Dopita (2002), because 
nitrogen is predominantly a primary nucleosyn- 
thesis clement in the range 12 + log(0/H) < 8.6 
(see also Fig. 3 in Kewley & Dopita (2002)) and 
the calibration of the (N/O) abundance ratio and 
the [Nn]/[0 n] line ratio is not sensitive to the the 
oxygen abundance under these circumstances. 

The results are summarized in Table 6 and 
Fig. 8 where we plot the three gradients as a func- 
tion of de-projected galactocentric radius. We find 
that inwards r~ 36" (~4.4 kpc) the O/H ratio 
follows a exponential profile with a slope of 0.140 
± 0.016 (dex/kpc) and 12 + log(O/H) r=0 ~ 8.9, 
similar to the normalized radial gradient found in 
other spiral disks. The outer abundance trend 
flattens out to an approximately constant value 
of 12 + log{O/H) r=0 ~ 8.27 (with a slight gradi- 
ent of 0.002 ± 0.019 (dex/kpc)) and could even re- 
verse (see Section 6 for a discussion on the possible 
causes for such flattening). The abundance gradi- 
ent derived, — 0.035dex/kpc (— O.OO^dex / arcsec) , 
is somewhat shallower than the one for the Milky 
Way (-0.08dex/kpc; Boissier & Prantzos 1999 
and references therein) . A similar trend was found 
from the analysis of the spectra of the concentric 
annuli. The analysis of these spectra was lim- 
ited to the innermost nine annuli because beyond 
this ring both the continuum and line emission 



become too faint to derive reliable emission-line 
fluxes. In Fig. 8 we plot the change in abundance 
for all annuli with a solid black line (bottom panel) 
where the uncertainties associated are represented 
by the grey-shaded area. The blue-dashed line 
represents the metallicity profile predicted by the 
best-fitting model of Boissier & Prantzos (1999, 
2000; hereafter BP2000 models). This profile has 
been shifted by 1 dex to match the metallicity scale 
derived from our spectroscopic data. Note that 
this kind of zero-point offsets in metallicity are not 
unexpected due to the significant uncertainties in 
the yields used in the disk evolution models (see 
discussion in Munoz-Mateos et al. 2011); besides, 
even the empirical oxygen abundances can be sub- 
ject to large systematic offsets (Moustakas et al. 
2010). The relative changes in metal abundances 
and therefore their radial profiles, on the other 
hand, are much more robust to these unknowns. 

5.4. Galaxy disk modeling 

In order to gain further insight into the evo- 
lution of NGC 5668, we have fitted its multi- 
wavelength surface brightness profiles with the 
BP2000 models. These models describe the chem- 
ical and spectro-photometric evolution of spiral 
disks as a function of only two variables: the di- 
mensionless spin parameter, A, and the circular 
velocity in the flat regime of the rotation curve, 
Vq- Within these models, galactic disks are simu- 
lated as a set concentric rings that evolve indepen- 
dently one from each other (for simplicity, radial 
mass or energy flows are not considered). The 
gas infall rate at each radius decreases exponen- 
tially with time, with a time-scale that depends 
on both the total mass of the galaxy and the lo- 
cal mass surface density at that radius. Once in 
the disk, gas is transformed into stars following 
a Kennicutt-Schmidt law multiplied by a dynami- 
cal term, which accounts for the periodic passage 
of spiral density waves. The mass distribution 
of each new generation of stars follows a Kroupa 
(2001) initial mass function (IMF). The finite life- 
times of stars of different masses is taken into ac- 
count; when they die, they inject metals into the 
ISM, thus affecting the metallicity of subsequent 
stars. The local metallicity at the time of forma- 
tion is taken into account when determining the 
lifetimes, yields, evolutionary tracks and spectra 
of each generation of stars. The model was first 
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calibrated against several observables in the Milky 
Way (see Boissier & Prantzos 1999) and then ex- 
tended to other galaxies with different values of A 
and Vq (Boissier & Prantzos 2000), using several 
scaling laws derived from the A-Cold Dark Matter 
framework of galaxy formation (Mo et al. 1998). 

For each pair of values of A and Vq, the model 
outputs radial profiles at different wavelengths, 
which can be then compared to the actual profiles 
of our galaxy. In order to probe the spatial loca- 
tion of stars of different ages —and therefore bet- 
ter constrain the model predictions— we measured 
surface brightness profiles at all GALEX FUV and 
NUV bands, the ugriz bands from SDSS and the 
3.6 and 4.5 /xm Spitzer bands. The radial vari- 
ation of internal extinction in the UV was esti- 
mated indirectly from the FUV— NUV color pro- 
files, as explained in section 5.1.1. The UV ex- 
tinction was then extrapolated to the optical and 
near-IR bands, assuming a Milky Way extinction 
curve, convolved with a sandwich model to ac- 
count for the relative geometry of stars and dust 
(see Munoz-Mateos et al. 2011 for details). 

We compared these extinction-free profiles with 
the model predictions for a grid of values of A and 
Vc , and used a \ 2 minimization algorithm to find 
the best fitting values. The results of the fit are 
shown in Fig. 9 and in Fig. 10 we present the 
residuals. Each panel corresponds to a different 
wavelength. The gray profiles are only corrected 
for foreground Milky Way extinction, whereas the 
black ones are corrected for internal extinction 
as well. The red and blue dashed lines bracket 
the region of the profile used for the fit. On 
one hand, starlight inside r ~ 20" is dominated 
by emission coming from the bulge and the oval. 
On the other hand, beyond r ~ 125" the signal- 
to-noise ratio becomes very low, and contamina- 
tion from background sources could be an issue. 
The model that best reproduces simultaneously all 
multi- wavelength profiles is shown with a red line. 
The shaded band encompasses all models whose 
X 2 is less than twice the x 2 °f the best model 
(Xmin -0.8) 14 . The resulting X 2 of the fitting 



Note that the scatter of the data-points around the model 
is not entirely due to noise, but also to real features of the 
galaxy such a smooth model cannot reproduce. This pre- 
cludes applying the classical statistical formula to translate 
the x 2 distribution into confidence intervals. The criterion 
of using x 2 < "^Xmin is j ust orientative, based on the re- 



model is obtained after a two-stage fitting pro- 
cedure because, as one would expect, this model 
does not reproduce the very small-scale variations 
of the surface brightness profiles of our disk. So 
in a first run of the fitting procedure we assume 
that the total uncertainty for each point is due 
to the quadratic sum of the zero-point and pho- 
tometric errors, plus an extra uncertainty in the 
model predictions of 10%. These values represent 
the initial guess for the intrinsic error of the model 
and in this stage the \ 2 is ~5. In the second step 
we calcute the rms of the best-fitting model with 
respect to the galaxy profiles and we pass these 
error profiles to the code as initial uncertainties to 
start the second run. In this way the new reduced 
X 2 values are close to unity. In principle, when 
we fit a set of data-points with a given model, 
it is often implicitly assumed that deviations be- 
tween the observed data-points and the model are 
due to the measured uncertainties of the former. 
However, in practice one also has to account for 
the fact that the models themselves are not a per- 
fect representation of nature and have their own 
"uncertainties" . The model yields smooth profiles 
that, by construction, cannot reproduce the fine 
structure of the actual profiles. In order to account 
for this, we performed for NGC 5668 this fit in a 
two-stage fashion, see also Sec. 5 of Munoz-Mateos 
et al. (2011) for a detailed discussion. When fitting 
the multi- wavelength profiles we should not ignore 
the internal degeneracy in the determination of the 
A and Vc values as these are not completely inde- 
pendent parameters. For this reason in Fig. 11 
we show the two-dimensional x 2 distribution ob- 
tained for NGC 5668 in the case of fitting all bands 
simultaneously and in the case of each band sep- 
arately. The corresponding best-fitting values are 
A = 0.053±g;gi| and V c = 167±g 2 kms" 1 . Where 
the errors are those expected when one is inter- 
ested in deriving each of the two quantities sepa- 
rately. 

In general, all bands are succesfully reproduced 
by the BP2000 models but the fit is not equally 
good at all wavelengths. For example, we can ap- 
preciate in the top pannels of Fig. 10 that in the 



suits shown in Fig. 9. Note also that we are not using the 
X 2 distributions in a strict statistical way, since many of 
the required mathematical conditions (gaussianity of er- 
rors, etc) are not met (see Munoz-Mateos et al. (2011) for 
a more extensive discussion). 
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cases of FUV, NUV and u-band data the quality 
of the fit is not so good because these bands are 
more sensitive to the recent variations in the star 
formation history and on the recipe used to calcu- 
late the exctinction (or even limitations intrinsic 
to the models in order to reproduce the UV part of 
the spectrum at certain metallicities, as explained 
in Fig.9 of Munoz-Mateos et al. (2011)). The de- 
viations in UV bands from the BP2000 models 
are compatible with the typical range of errors, 
as pointed out in Munoz-Mateos et al. (2011). 

While the BP2000 models provide a good fit to 
the overall shape of the surface brightness profiles 
of NGC 5668 (see Fig. 9), the analysis of the devi- 
ations from these otherwise idealized models could 
give important clues on the details of the star 
formation history of this galaxy (e.g. gas and/or 
stellar radial transfers). In this regard, Fig. 12 
shows that while the best-fitting model predicts a 
systematic bluing in the colors towards the outer 
parts the disks, the measured colors flatten or even 
get redder beyond ^30-40". This radius interest- 
ingly coincides with the position where the slope 
of the metallicity gradient changes sign. We also 
find small offsets (<0.2 mag) between the observed 
and predicted colors. This is partially due to un- 
certainties in the model predictions for some of 
these colors. 

Note that the model has to reproduce at the 
same time the shape of the observed surface 
brightness profiles in all bands simultaneously and 
in addition it has to reproduce the "average" level 
in surface brightness. This kind of fitting intro- 
duces some limits in the range of the model pa- 
rameters in a way that the colors of these models 
are located in a relatively narrow range as a con- 
sequence of the high degeneration within those 
model (surface brightness) profiles that match the 
data at all wavelengths simultaneously . 

Small variations of a few tenths of a magnitude 
will not have a significant impact on the overall 
shape of the surface brightness profiles, but will 



One could decide to fit only the color profiles instead. How- 
ever, that would lead to solutions that might not fit at all 
the overall shape of the surface brightness profiles and that 
would be very sensitive to uncertainties in the model pre- 
dictions regarding specific colors (associated to limitations 
in the yields, stellar libraries at non-solar abundances, un- 
certainties in the luminosity at near-infrared and UV wave- 
lengths; sec Maraston et al. (2006) and Munoz-Mateos et al. 
(2011), respectively). 



reveal themselves more clearly in the color pro- 
files (see Munoz-Mateos et al. 2011). On the other 
hand, the SFH of the different annuli in the model 
are tied one with each other (due to the analyt- 
ical way in which we implement radial changes) 
but the actual SFH within the galaxy will exhibit 
more complex radial variations from one ring to 
the next one. 

6. Discussion 

The analysis presented above shows that most 
of the properties of the disk of NGC 5668 can 
be well reproduced by the chemical and spectro- 
photometric BP2000 models. This is true both for 
the global metallicity gradient as well for the shape 
of the galaxy surface brightness profiles. This 
is also true for the rotation curve (see Fig. 13). 
However, there are observables (chemical abun- 
dances and colors) and, especially, radial inter- 
vals among them that are not so well fitted by 
these simple models. The objective must then be 
to establish the mechanism(s) that drive the de- 
partures of these observables and the associated 
star-formation and chemical histories in NGC 5668 
from those predicted by the inside — out scenario 
of disk formation on which the BP2000 models are 
based on. In particular, we ought to explain the 
bimodality (Chiappini et al. 2001) of the abun- 
dance gradient (and color profiles) measured in 
NGC 5668: flatter in the outermost regions and 
steeper in the inner ones. 

6.1. Bimodal distribution mechanisms 

The re-distribution of stars, gas and dust in 
galaxies can have important consequences in the 
metallicity distribution and, in particular, in de- 
viations of the chemical and photometric prop- 
erties of disks with respect to the predictions of 
static disk-evolution models. In this regard, radial 
changes (up— or down— bending) in the slope of 
the abundance gradient of nearby spiral galaxies 
(Vilchez & Esteban 1996) are detected and are 
often related to gaseous and stellar radial mix- 
ing processes (Spitzer & Schwarzschild 1953; Bar- 
banis & Woltjer 1967; Shaver et al. 1983; Fuchs 
2001; Sellwood & Binney 2002; Roskar et al. 2008; 
Haywood 2008; Bresolin et al. 2009b; Vlajic et al. 
2009). 

The existence of this kind of bimodal gradients 
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Fig. 9. — Best-fitting model of NGC 5668 using the Kroupa (2001) IMF. The gray profiles are only corrected for foreground 
Milky Way extinction, whereas the black ones are also corrected for internal attenuation. The latter profiles are the ones used 
to constrain the disk-evolution models in the radial range spanned between the red and blue vertical lines (in order to exclude 
the bulge and the low S/N outer parts). The best-fitting model is shown as a red solid line, and the red shaded band comprises 
all models with x 2 < 2x^ in . The radius along the semimajor axis is shown both in kpc (bottom x axis) and arcseconds (top x 
axis). 
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Fig. 10. — Residuals of the best-fitting BP2000 model of NGC 5668 using the Kroupa (2001) IMF. The red line represents 
the best- fitting model and the red shaded area represents all models with x 2 5: 2x^ in . The black points are the residuals of 
the black profile corrected for internal and foreground Milky Way attenuation. The radius along the semimajor axis is shown 
both in kpc (bottom x axis) and arcseconds (top x axis). 
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Fig. 11. — Two-dimensional x 2 distribution of the disk-evolution model Biting procedure. In the top-left panel we plot x 2 as 
function of the spin (A) for different circular velocities (Vc). The best-fitting value and its estimated uncertainty are marked 
with a circle with error bars. The middle-left plot is analogous and shows x 2 as a function of Vc for different spins. These are 
projections of the surface shown in the bottom-left panel, where the red dot marks the best model. While we constrain the 
model by fitting all bands at the same time, we keep track of the individual x 2 distributions at each wavelength. These are 
shown in the small panels to the right (the red point has been replicated in all of them for clarity) . The missing small panels 
would correspond to J, H, and K 2MASS images that we decided to exclude from this analysis as they are not deep enough to 
provide additional information. For more details see Munoz-Mateos et al. (2011). 



25 



r (arcsec) 

50 100 150 200 




5 10 15 20 25 30 

r (kpc) 



Fig. 12. — Observed color profiles of NGC 5668 (points) corrected for internal extinction as in the black points of Figure 9. 
These observed colors are compared with the ones predicted by the best- fitting model to light profiles (lines and bands). The red 
and blue vertical dashed lines bracket the spatial range used during the fit. Note that since the fit is performed simultaneously 
at all radii and all wavelengths from the FUV to 4.5 /im, departures up to a few tenths of a magnitude are expected for a 
particular color and radius. 
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Fig. 13. — Left: comparison between the observed and predicted rotation curve of NGC 5668. The black solid line shows the 
rotation curve corresponding to the model that best fits the multi-wavelength profiles, and the gray band shows the associated 
uncertainty. The empirical data come from the Fabry-Perot Ha studies of Jimenez- Vicente & Battaner (2000) and Epinat et 
al. (2008). The rotation curve of the former is plotted with a dashed line, the dotted ones showing the corresponding errors. 
The values of Epinat et al. (2008) are marked with solid points (receding side of the disk) and open ones (approaching side). 
Since NGC 5668 is almost face-on, small variations in the adopted inclination angle could rescale the empirical rotation curves. 
The vertical tickmark near the end of the curve shows the variations induced by changes of±l° in the inclination angle. Right: 
Ho velocity map of the ionized gas derived from the Ha emission line (Ha is represented by isocontours). North is up and East 
is to the left. One pointing was masked out because the derived velocity values were not reliable. 



are often proposed to have been induced by the 
presence of a bar-like potential (Diaz et al. 1990; 
Vila-Costas & Edmunds 1992; Edmunds & Roy 
1993; Zaritsky et al. 1994). Indeed, in the case of 
strong-barred galaxies a change from a shallow to 
a steep metallicity profile is commonly observed 
and explained in the context of radial mixing of 
the gas induced by the bar (Martin & Roy 1995; 
Roy & Walsh 1997). On the other hand, bars can 
also lead to a steepening of the metallicity profile 
towards the inner disk if the gas transfer results 
in significant star formation in situ (Friedli et al. 
1994; Martin & Roy 1995; Roy et al. 1996; Roy & 
Walsh 1997). 

This latter behavior has been found in dynam- 
ical simulations of the formation of bars (Friedli 
et al. 1994; Friedli & Benz 1995). In these sim- 
ulations the presence of a steep-shallow break in 
the metallicity profile is the result of an intense 
chemical enrichment by star formation in the bar 
combined with the dilution effect of the outward 
flow beyond the break. According to these models 
the presence of such a break indicates that the bar 
has recently formed, i.e. in the last Gyr. Observa- 



tionally, the metallicity distribution of a number 
of galaxies hosting young bars have been success- 
fully reproduced by this or a similar scenario (Roy 
& Walsh 1997; Considere et al. 2000). 

However, other mechanisms have been also pro- 
posed to explain bimodal metallicity distributions. 
Thus, in the context of the spiral density- wave 
theory, star formation is expected to be propor- 
tional to 51 — n, p (Oort & Peixoto 1974). In Jensen 
et al. (1976) this prediction was successfully tested 
against the metallicity gradients measured in a 
number of grand-design spirals. Should this sce- 
nario be valid one would expect to find a mini- 
mum in the chemical abundance at the position of 
the corotation radius (see also the recent works by 
Acharova et al. 2005 and Scarano 2010), especially 
in the case of galaxies where the spiral structure 
is long-lived and quasi-stationary (McCall 1986; 
Acharova et al. 2005). This is to the fact that 
f2 — flp is null (by definition) at corotation and, 
therefore, star formation should less efficient re- 
sulting in a lesser degree of metal enrichment. In 
this scenario a maximum in the colors should be 
also observed at the approximate position of the 
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corotation radius as the optical-near-infrared col- 
ors would be dominated by those of the underlying 
stellar population since the current day SFR is ex- 
pected to be low. 

Finally, a flattening of the outer color and 
metallicity profiles is believed to be also caused by 
the increasing contribution of old migrated stars 
to the stellar content of disks (Zaritsky et al. 1994; 
Binncy 2001; Scllwood & Binney 2002; Roskar 
et al. 2008 and more recently Minchev et al. 2011; 
Roskar et al. 2011). In this case, as we move 
further out into the outer disk we find progres- 
sively older stars, which have been migrating for 
a longer period of time, leading to a positive age 
and color gradient in these regions. With regard to 
the metal abundance and according to the simula- 
tions of Roskar et al. (2008) radial migration leads 
to the mixing of the old stellar populations, which 
results in flatter gradients at early times and in 
the very outer regions of present-day disks as these 
regions are primarly populated by migrated stars. 
In the few observational and theoretical studies on 
the luminosity- weighted age profiles of the disks a 
similar flattening or even an up-bending is found 
(Bakos et al. 2008; Sanchez-Blazquez et al. 2009; 
Yoachim et al. 2010). While some authors favor 
stellar migration as the main driver of the shaping 
of these age profiles (Yoachim ct al. 2010) some 
others indicate that this might due to a decrease 
in the star formation in the external parts of the 
disk with time caused by a reduction of the vol- 
ume density of the gas in these regions (Sanchez- 
Blazquez et al. 2009). 

The different scenarios proposed above to ex- 
plain a steep-shallow break in the metallicity gra- 
dient of galaxies, namely a young bar, reduced star 
formation at the corotation radius, stellar migra- 
tion or evolution of the star-formation threshold 
with rcdshift have different imprints not only on 
the specific shape of the metallicity gradient but 
also in other properties such as the colors. In 
the case of NGC 5668 the change in the shape 
of the metallicity gradient takes place at a ra- 
dius of ^36"(4.4kpc or 2.8 disk-scale lengths), 
well within the region where in-situ star formation 
takes place in the disk, at an approximate surface 
brightness of ^22magarcsec~ 2 (3.6 (im Spitzer 
band). At this radius the contribution of migrated 
stars to either the colors, luminosity-weighted age, 
or chemical abundance is expected to be negligible 



compared with that from stars formed in-situ (see 
Bakos et al. 2008, Vlajic et al. 2009). Again, as 
this brakes takes place well within the star-forming 
disk of NGC 5668 changes in the color and metal- 
licity gradients associated with a possible evolu- 
tion of the star-formation threshold with rcdshift 
is highly unlikely. With respect to the possibility 
that the minimum in the metal abundance profile 
of NGC 5668 could be due to the presence of the 
corotation radius at this break two things can be 
said. Firstly, NGC 5668 is a rather flocculent spi- 
ral 16 , where spiral density waves are expected to 
be weak or absent and, consequently, the effects 
of the spiral arms on the radial distribution of the 
star formation should be minimal (if any; McCall 
1986). Secondly, the color profiles in (u — g) and 
(g — r) show a minimum at the position of the 
metallicity break, which is the opposite to what 
we would expect if that position corresponds to 
the corotation radius. 

6.2. Bar formation in NGC 5668 

The only scenario that remains to be analyzed 
in detail is the possibility that the deviations of the 
metallicity and color profiles from those predicted 
by the best-fitting BP2000 model of NGC 5668 are 
due to the presence of a nascent bar and the ef- 
fects associated to it. It is well known that one 
of the most important drivers in the evolution 
of galaxies are bars (or non axisymmetric central 
light distributions or ovals, Athanassoula 1994). 
Bars exist in about two-thirs of disk galaxies (Scll- 
wood & Wilkinson 1993) and isolated galaxies are 
known to be able to develop a barred morphology 
spontaneously from internal instabilities caused 
by cooling processes (Miller & Prendergast 1968; 
Kalnajs 1978; Binney & Tremaine 1987; Sellwood 
& Wilkinson 1993). These gravitational instabili- 
ties are commonly the result of enhanced gas ac- 
cretion in the disk (Lindblad et al. 1996). Some of 
the observational properties of NGC 5668 already 
reveal a potential active gas-accretion phase in this 
galaxy, such as the presence of both HVCs and 
HRVCs (Jimenez- Vicente & Battaner 2000 and 
references therein) and the high total star forma- 
tion rate in this object. 



The lack of strong spiral arms in NGC 5668 is evident even 
after the examination of the IRAC near-infrared images of 
the object (see Scarano 2010 and references therein). 
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In this regard, another relevant parameter is the 
HI content of NGC 5668 compared with objects of 
similar total mass and morphological type. In or- 
der to find out whether the HI content NGC 5668 
is particularly high or low, we resort to the so- 
called HI deficiency parameter. This quantity, de- 
fined by Haynes & Giovanelli (1984), compares 
—in logarithmic scale— the observed HI mass of 
a given galaxy and the typical HI mass of isolated 
galaxies of similar morphological types and linear 
sizes. The difference between both values is per- 
formed in such a way that positive HI deficiencies 
correspond to galaxies with less gas than similar 
field galaxies, and viceversa. Following the pre- 
scriptions of Haynes & Giovanelli (1984), we find 
that an Sd galaxy with the same optical diameter 
as NGC 5668 is expected to have log(Mni) = 9.66. 
Solanes et al. (1996) extended the work of Haynes 
& Giovanelli (1984) to a larger sample of galax- 
ies. While they only include Sa-Sc galaxies in 
their sample, we can safely apply the fitting co- 
efficients of Sc's to our Sd spiral (Solanes et al. 
2001). By doing so we obtain a reference HI mass 
of log(Mni) = 9.64, completely consistent with 
the previous estimation. According to Schulman 
et al. (1996), the actual HI mass of NGC 5668 is 
log(Mni) = 10, which implies a negative HI de- 
ficiency of ~ —0.35. Considering that the usual 
la scatter of the HI mass for a given type and 
size is 0.2—0.3, dex (including all types of galax- 
ies), we conclude that our galaxy is roughly l-2cr 
gas-richer than its normal counterparts. While 
this difference might not be significant it shows 
that if anything, NGC 5668 has a larger HI con- 
tent than the average of the spiral galaxy popula- 
tion of its type. Schulman et al. (1994) found in 
NGC 5668 a weak bar or oval inner structure of 
12" size visible in both optical and near-infrared 
images. Athanassoula (1994) found that in spiral 
galaxies having oval structures or weak bars star 
formation can occur prolifically along them, espe- 
cially in the central parts and at the two ends of 
the bar region. This would result in bluer colors in 
these regions. The presence of a local minimum in 
color at a galactocentric distance of ~30-40 arcsec 
seems to favor this scenario for NGC 5668. Note, 
however, that the oval seen in its images is signifi- 
cantly smaller than the radius where the minimum 
in color and metallicity is found. One possible ex- 
planation for this difference might come from the 



fact that, while the light from the superposition of 
xi orbits that shape the oval is only clearly seen 
in the central 12", the instabilities associated to it 
could take place further out, where the superpo- 
sition of these orbits is not yet significant enough 
for being detected via photometry only (Debat- 
tista et al. 2006). 

Another possible explanation for the existence 
of the oval distortion in NGC 5668 is the possible 
interaction of NGC 5668 with UGC9380 Schulman 
et al. (1996). This is a dwarf galaxy at a relative 
distance of 200 kpc to the southeast of NGC 5668, 
with a systemic velocity of 1690 km sec -1 and a 
relative velocity of ^108 kmsec -1 . A recent or on- 
going tidal interaction between these two galaxies 
could have produced the high velocity features in 
NGC 5668 and triggered the formation of both the 
oval and the incipient bar, as reported by N-body 
simulations of minor mergers Eliche-Moral et al. 
(2006) and Eliche-Moral et al. (2011). 

The fact that these blue colors do not extend 
all the way from this galactocentric distance to the 
center of the galaxy is likely due to a combination 
of different effects. Firstly, the secular inside-out 
model for the evolution of the disk of NGC 5668 
predicts a reddening of ^0.2 mag (0.15 mag) in u — 
g (g — r) from 36 to 6" (see Fig. 12). Moreover, the 
reddening in the (FUV — NUV) color indicates a 
change of ^0.09 mag in the B — V color purely due 
to dust attenuation [E(B — V)] within this same 
radial range. Finally, the light contribution of the 
bulge in the center, despite being small, can also 
lead to relatively red colors. These three effects arc 
supcrpimposed on the change in colors induced by 
star formation in the bar region and whose effects 
on the colors appear to be noticeable only at both 
its ends. 

Regarding the age of this young bar, in the case 
of the barred galaxy NGC 3359, Martin & Roy 
(1995) calculated the age of bar by using the equa- 
tion of turbulent transport in a shear flow (Roy & 
Kunth 1995). Below we repeat this exercise for 
the case of NGC 5668. If we consider a cartesian 
coordinate system centered in the galaxy center 
and we take x 2 as the radial direction of the local 
circular orbit (we do not take in account perpen- 
dicular effects), I ^300 pc as the mean free path 
of clouds between collisions (Roberts & Hausman 
1984), we can then calculate the time for the gas 
to diffuse in a length Ax2 in the radial direction 
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as: 



agreement with the predictions of Calzetti 
(2001). 



In our case Ax2 = 4.4 kpc, wich coincides with 
t break is the radius at which we detect the flat- 
tening of the metallicity gradient and the mini- 
mum in the color profiles. From the radial- velocity 
map kindly provided by J. Jimenez- Vicente (priv. 
comm. ; see also Jimenez- Vicente & Battaner 2000) 
we determine the maximum of the non-circular 
motions to be v—10 kms^ 1 . This yields t X2 ~ 10 9 
yr, which can be considered as a upper limit to 
the bar age since the radial transport induced by 
a bar is a stationary flow. The lower limit is given 
by simple division of Aa^/v. The best estimate 
for the age of the young bar in NGC 5668 is then 
10 8 < t X2 < 10 9 yr. 

Summarizing, while the overall observational 
properties of NGC 5668 are well fitted by the 
inside-out scenario of disk formation, the devia- 
tions in the color and metallicity profiles are best 
interpreted in the context of the presence of a 
nascent bar where significant in-situ star forma- 
tion is (or have been recently) taking place. The 
formation of the bar is believed to be due to in- 
stabilities in the gas-rich inner disk of NGC 5668 
possibly helped in the interaction with a compan- 
ion. This scenario is compatible with the rela- 
tively large HI content and star formation rate of 
NGC 5668 and with the presence of HVCs and 
HRVCs in its velocity field, evidence of significant 
non-rotational gas motions in the disk. 

7. Conclusions 

In this paper we have carried out an exten- 
sive and detailed study of the chemical and pho- 
tometric properties of the nearby spiral galaxy 
NGC 5668. This detailed study has been possi- 
ble thanks to the combined use of Integral Field 
Spectroscopy (IFS) in the optical and panchro- 
matic broad-band imaging of the entire system. 
The main conclusions from this work are the fol- 
lowing: 

1. Dust- attenuation profiles have been ob- 
tained using a number of methods. The 
mean continuum attenuation is Ay ~0.4 mag. 
We find a significantly larger ionized-gas at- 
tenuation than that of the continuum in 



2. We have derived dust- attenuation-corrected 
emission-line fluxes for a total of 62 indi- 
vidual H II complexes and 18 concentric 
annuli centered on the position of the nu- 
cleus of NGC 5668. Based on the strong- 
line method of Kewley & Dopita (200X) 
oxygen-abundance measurements have been 
obtained. We find a bimodal radial distribu- 
tion in metallicity with a steep negative gra- 
dient with slope -0.14±0.09(dex/kpc) within 
a galactocentric distance of 36"and a shal- 
lower (or even positive) metallicity gradient 
beyond that point, 0.002±0.002 (dcx/kpc). 

3. Surface brightness and color profiles from 
the UV to the near-infrared have been also 
obtained. These profiles nicely match the 
predictions of the chemo-spectrophotometrical 
models for the evolution of galaxy disks 
(BP2000 models) for a circular velocity of 
Wcirc = 167 km s^ 1 and a spin parameter of 
A=0.053. This best-fitting model also agrees 
with the overall shape of the galaxy metal- 
licity gradient and rotation curve although 
it cannot reproduce the steep-shallow metal- 
licity break neither the moderate bluing in 
some of the color profiles (u—g, g — r) around 
the position of the metallicity break. 

4. Out of the different mechanisms proposed 
in the literature to explain the change in 
the slope of the metallicity profile in spiral 
disks (and the color profile), only the pres- 
ence of a bar in its formative stages agrees 
well with the other observational properties 
of NGC 5668: the position of the metallicity 
break in disk-scale lengths, the bluing in the 
colors at that position, the large HI content 
and star formation rate of the galaxy, the 
presence of HVCs and HRVCs in the galaxy 
velocity field and even the detection of an 
oval in the central region of the galaxy. 
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Tabic 1: Global properties of NGC 5668. 



Name 



NGC 5668 



Source 



Morphological type 

RA(2000) 

Dec(2000) 

B 

R 

L B 

Lfir 
D 25 

Heliocentric systemic velocity 

Distance 

Inclination 

PA 

Total dynamical mass 
E(B-V) 



SA(s)d 
Sc(s)II-III 
14h 33m 24.3s 
+04° 27' 02" 
12.13 ± 0.03 mag 
11.28 ± 0.01 mag 
2.7 ± 0.6 10 9 L Q Q 
1.0 ±0.3 x 10 8 L Q a 
(5.8 ± 1.2) x 10 9 L Q a 
3.3' ± 0.2' 
1582 ± 5 kms- 1 
24.8 ± 1.7 Mpc 
18° 
145° 
5.7 x 10 10 Mq 
0.037 mag 



de Vaucouleurs et al. (1991) 
Sandage & Tammann (1987) 
Dressel & Condon (1976) 
Dressel & Condon (1976) 
Schulman (1997) 
Schulman (1997) 
Schulman (1997) 
Schulman (1997) 
Schulman (1997) 
Jimenez & Battaner (2000) 
Schulman (1997) 

Fixen (1996) 
Schulman (1997) 
Schulman (1997) 
Schulman (1997) 
NED 



a Corrected for Galactic and internal extinction. 

6 Uncorrected for Galactic and internal extinction L Qi b = 5 x 10 32 ergs -1 



Table 2: NGC 5668 Observational Log. 



Pointing^ 


Obs. Date 


Offsets 


Exposure Time 


Air-mass 


Seeing 




(UT) 


(arcsec) 


(sec) 




(arcsec) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1 


2007-06-22 


(0,0) 


3x1000 


1.21 


1.1 


2 


2007-06-22 


(0,60) 


4x1000 


1.38 


1.1 


3 


2007-06-22/23 


(52,60) 


3x1000 


1.75 


1.1 


7 


2007-06-23 


(52,-60) 


3x1000 


1.21 


1.2 


5 


2007-06-23 


(0,-60) 


3x1000 


1.36 


1.2 


4 


2007-06-23/24 


(-52,60) 


3x1000 


1.74 


1.0 



NOTE: (1): NGC 5668 pointings. (2): Observations date. (3): Pointings offsets; the original pointing of 
the telescope is at RA(2000) 14 h 33 m 24.3 s and Dec(2000) +04°33'24.3". (4): Exposure time in seconds and 
numbers of images obtained for each pointing. (5)-(6): Mean air-mass and seeing values, 
t Pointing number 6 in the default mapping strategy used with PPAK at the CAHA observatory has offsets 
(-52", -60") relative to the center of the mosaic and was not observed in the case of NGC 5668. 



36 



Table 3 
Spectra emission lines. 



ID 


F[OJ/]/FHa 


F[OIII]/Fua 


F[AT//]/F H a 


F[S//a]/FHa F[5// b ]/F H Q 


FHa 


Ay 


(1) 


(2) 


(3) 


(4) 


(5) (6) 


(7) 


(8) 



Rings 



1 
2 
3 
4 
5 
G 
7 
8 
9 

10 
11 
12 
13 
14 
15 



2 


68±0 


43 








0.32±0.05 


0.28±0.05 


0.27±0.05 


18 


20±2 


51 


1.57 


1 


92±0 


31 





15±0 


04 


0.31±0.05 


0.26±0.04 


0.24±0.04 


19 


40±2 


41 


1.52 


1 


77±0 


29 





21±0 


05 


0.26±0.04 


0.25±0.04 


0.20±0.04 


13 


94±1 


77 


1.23 


1 


13±0 


21 





21±0 


05 


0.29±0.06 


0.25±0.06 


0.23±0.05 


4 


96±0 


71 


0.35 


1 


70±0 


32 





28±0 


06 


0.26±0.05 


0.25±0.05 


0.19±0.05 


9 


36±1 


41 


1.00 


1 


83±0 


24 





36±0 


05 


0.21±0.03 


0.23±0.03 


0.18±0.03 


10 


61±1 


18 


1.33 


1 


30±0 


18 





40±0 


06 


0.17±0.03 


0.21±0.03 


0.17±0.03 


7 


12±0 


78 


0.88 


2 


29±0 


41 





47±0 


09 


0.20±0.05 


0.24±0.05 


0.25±0.05 


7 


86±1 


29 


1.64 


2 


04±0 


48 





61±0 


14 


0.15±0.06 


0.23±0.07 


0.19±0.06 


4 


36±0 


93 


1.33 


1 


81±0 


58 





35±0 


13 




0.21±0.09 


0.21±0.09 


2 


80±0 


80 


1.19 


1 


99±0 


87 





64±0 


28 








1 


80±0 


75 


1.10 


1 


50±0 


43 





66±0 


18 








2 


03±0 


53 


0.76 



HII Regions 



1 


1 


14±0 


14 





64±0 


04 





09±0 


02 





28±0 


03 





19±0 


03 


2.23±0 


09 


0.65 


2 


1 


19±0 


37 





48±0 


11 





20±0 


09 














0.50±0 


05 


0.12 


3 


1 


55±0 


23 





49±0 


04 





15±0 


02 





24±0 


02 





13±0 


01 


11.36±0 


28 


2.02 


4 


1 


52±0 


09 





24±0 


02 





16±0 


01 





21±0 


01 





17±0 


01 


12.27±0 


16 


1.36 


5 


1 


33±0 


06 





37±0 


01 





17±0 


01 





22±0 


01 





15±0 


01 


13.32±0 


14 


1.18 


6 

















01±0 


01 





28±0 


05 





25±0 


05 


5.39±0 


43 


2.46 


7 


1 


83±0 


42 





44±0 


08 





02±0 


01 





28±0 


05 





13±0 


03 


3.68±0 


21 


1.64 


8 


1 


44±0 


12 





29±0 


04 





09±0 


03 





18±0 


03 





21±0 


04 


11.39±0 


50 


1.31 


9 


1 


57±0 


11 





46±0 


03 





12±0 


02 





12±0 


02 





10±0 


02 


32.00±0 


91 


2.06 


10 


1 


72±0 


07 





41±0 


01 





17±0 


01 





17±0 


01 





13±0 


01 


12.39±0 


16 


1.22 


11 


2 


00±0 


22 





62±0 


05 





14±0 


02 





26±0 


02 





20±0 


02 


4.75±0 


18 


1.30 


12 


2 


05±0 


26 





67±0 


04 





11±0 


02 





24±0 


02 





14±0 


02 


10.02±0 


30 


1.98 


13 


2 


31±0 


23 





34±0 


06 





18±0 


05 





26±0 


06 





14±0 


05 


4.27±0 


25 


0.39 


14 


1 


76±0 


80 





32±0 


10 





20±0 


05 





16±0 


06 








23.48±1 


85 


2.60 


15 





48±0 


08 





57±0 


05 





10±0 


04 





15±0 


05 





12±0 


04 


6.27±0 


25 


0.73 


16 


1 


18±0 


09 





47±0 


02 





15±0 


01 





17±0 


01 





11±0 


01 


15.94±0 


23 


1.61 


17 


3 


93±0 


76 





96±0 


16 





18±0 


07 





22±0 


06 





22±0 


09 


19.32±1 


93 


2.89 


18 


1 


26±0 


37 











22±0 


03 





29±0 


03 





13±0 


03 


6.84±0 


30 


1.93 


19 


1 


48±0 


17 





29±0 


04 





26±0 


03 





27±0 


03 





18±0 


03 


11.43±0 


50 


0.90 


20 





51±0 


12 





41±0 


04 





25±0 


03 





28±0 


04 





17±0 


03 


1.375±0 


05 


0.12 


21 


1 


85±0 


14 





29±0 


03 





22±0 


02 





15±0 


02 





11±0 


02 


16.73±0 


59 


0.89 


22 


2 


40±0 


62 





26±0 


10 





18±0 


05 





35±0 


05 





24±0 


04 


4.30±0 


34 


1.87 
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Table 3 — Continued 



ID F[ jj]/F Hq F[ ///]/F H a F[jvjj]/Fh q F[s /Jo ]/F H a F[s/ /b ]/F HQ F Hq Ay 

(1) (2) (3) (4) (5) (6) (7) (8) 



23 





.93±0.04 





.97±0 


.01 





.09±0 


.01 





12±0 


.01 


0. 


10±0 


.01 


11.94±0 


.12 


0.92 


24 


1 


.69±0.15 





.26±0. 


.03 





,23±0 


.02 


0. 


22±0 


.02 


0. 


16±0 


.02 


7.36±0 


.18 


1.31 


25 


1 


,57±0.27 





,71±0. 


.07 





.24±0. 


.04 


0. 


14±0 


.03 


0. 


10±0 


.03 


2.41±0 


.12 


1.06 


26 


1 


,41±0.51 





14±0 


.05 





,31±0. 


.07 


0. 


33±0 


.07 





21±0 


.05 


9. Mil. 


.68 


0.88 


27 






2 


19±0 


.48 








0. 


15±0 


.07 








31.30±4. 


.64 


4.59 


28 


1 


.74±0.47 





17±0. 


.05 





30±0 


.05 


0. 


26±0 


.04 


0. 


,19±0 


.04 


20.57±2. 


.59 


1.67 


29 


1 


.28±0.20 





,21±0 


.03 





.25±0. 


.03 


0. 


24±0 


.03 


0. 


19±0 


.03 


14.89±0. 


.93 


1.08 


30 


1 


,34±0.24 





,24±0. 


.05 





,25±0. 


.04 


0. 


,25±0 


.04 


0. 


16±0 


.03 


11.61±0 


.71 


1.15 


31 


1 


.30±0.08 





.38±0. 


.02 





,18±0. 


.01 





18±0 


.01 





11±0 


.01 


9.23±0 


.12 


1.14 


32 


1 


13±0.09 





.33±0 


.02 





,18±0. 


,01 





15±0 


.01 


0. 


11±0 


.01 


29.66±0 


.59 


1.18 


33 


1 


.43±0.13 





.20±0. 


.02 





.24±0. 


.02 





19±0 


.02 





12±0 


.01 


30.84±1. 


.21 


1.30 


34 


2 


,30±0.59 





.41±0 


.10 





,14±0. 


.04 





,24±0 


.03 


0. 


13±0 


.03 


6.59±0 


,41 


2.20 


35 


3 


,42±0.38 





.39±0 


.06 





.35±0. 


.04 





35±0 


.04 





.21±0 


.04 


50.93±3. 


.45 


1.50 


36 


2 


,32±0.15 





.66±0 


.03 





,14±0 


.01 





19±0 


.01 


0. 


12±0 


.01 


17.41±0. 


.34 


1.96 


37 


1 


,15±0.05 





.36±0 


.01 





,18±0 


.01 


0. 


14±0 


.01 


0. 


11±0 


.01 


57.29±1. 


.04 


1.00 


38 


1 


.55±0.23 





15±0 


.04 





.26±0. 


.03 


0. 


24±0 


.03 


0. 


15±0 


.03 


13.94±0. 


.79 


1.25 


39 


2 


,68±0.92 





96±0 


.14 




















7.80±0 


.59 


2.74 


40 


3 


,32±1.16 





,51±0. 


.11 





15±0 


.03 


0. 


26±0 


.03 


0. 


,09±0 


.02 


25.61±1. 


.71 


3.63 


41 


1 


,06±0.05 





.32±0. 


.01 





,16±0. 


.01 





19±0 


.01 


0. 


15±0 


.01 


14.14±0. 


.16 


1.06 


42 


2 


,21±0.45 





.27±0. 


.07 





.24±0. 


.05 


0. 


.23±0 


.04 


0. 


.20±0 


.04 


7.85±1. 


.30 


1.84 


43 





,76±0.16 





,20±0. 


.03 


0. 


,22±0. 


.03 





23±0 


.03 


0. 


18±0 


.03 


18.11±1. 


.75 


1.22 


44 


1 


,78±0.19 





,15±0 


.03 





,32±0 


.03 


0. 


,23±0 


.02 





,24±0 


.02 


28.78±2. 


.28 


1.47 


45 


2 


,99±0.41 





12±0 


.04 





,29±0. 


.04 


0. 


.27±0 


.04 





,21±0 


.03 


25.34±3. 


.02 


1.75 


46 


1 


,33±0.25 





15±0. 


.05 





.25±0. 


.05 


0. 


26±0 


.05 


0. 


,21±0 


.04 


11.02±1. 


.50 


1.06 


47 


1 


,45±0.20 





16±0 


.04 





,25±0. 


.04 


0. 


26±0 


.04 


0. 


,21±0 


.04 


6.23±0. 


.71 


0.34 


48 


3 


,60±0.75 





.26±0. 


.08 





,26±0 


.06 





29±0 


.07 





18±0 


.05 


15.20±2. 


.89 


1.71 


49 


1 


,64±0.05 





,75±0. 


.01 





10±0 


.01 





11±0 


.01 


0, 


10±0 


.01 


16.48±0. 


.16 


1.09 


50 


1 


,22±0.19 





.20±0. 


.04 





.25±0. 


.04 


0. 


25±0 


.04 


0, 


21±0 


.04 


9.11±1. 


.05 


0.72 


51 


1 


,56±0.11 





.34±0. 


.03 





.21±0. 


.02 


0. 


20±0 


.02 





17±0 


.02 


22.45±1. 


.29 


0.93 


52 


1 


.14±0.11 





.31±0. 


.03 





,22±0. 


.02 





,22±0 


.02 


0. 


19±0 


.02 


21.20±1 


.41 


1.01 


53 


1 


,52±0.05 





.40±0. 


.01 





18±0. 


.01 





19±0 


.01 





16±0 


.01 


29.32±0. 


.54 


1.02 


54 


1 


,65±0.06 





.66±0. 


.02 





12±0. 


.01 





15±0 


.01 





12±0 


.01 


10.30±0. 


.12 


0.94 


55 


1 


,72±0.08 





,63±0. 


.02 





12±0. 


.01 


0. 


,22±0 


.01 





16±0 


.01 


15.11±0. 


.18 


1.33 


56 





.81±0.05 





18±0. 


.02 





,11±0. 


.02 


0. 


,24±0 


.02 


0. 


16±0 


.02 


3.30±0. 


.16 


0.00 


57 


11.5±2.4 





,78±0. 


.16 





18±0. 


.04 





18±0 


.04 


0. 


16±0 


.04 


153.29±30 


.27 


3.93 


58 


1.83±0.08 





.70±0. 


.02 





,11±0. 


.02 





15±0 


.02 


0. 


13±0 


.03 


6.07±0. 


.16 


0.86 


59 


3.22±0.25 





,45±0. 


.04 





,13±0 


.04 


0. 


19±0 


.03 


0. 


27±0 


.04 


6.09±0. 


.27 


1.51 


60 


2.26±0.14 





36±0 


.03 





,10±0. 


.03 


0. 


17±0 


.03 


0. 


16±0 


.04 


4.80±0. 


.18 


1.06 


61 


7.2±1.3 





.78±0. 


.11 








0. 


17±0 


.05 


0. 


,21±0 


.07 


20.34±1. 


.59 


3.57 
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Table 3 — Continued 



ID F[0//]/FHa F[o///]/FHa F[jv//]/Fh q F[SJIo] /Fho F[5// b ]/FHa Fna Ay 

(1) (2) (3) (4) (5) (6) (7) (8) 



62 3.24±0.61 ••• 0.28±0.05 0.25±0.05 0.21±0.04 34.30±5.43 2.30 

Total* 1.31±0.17 0.35±0.05 0.16±0.04 0.21±0.03 0.15±0.03 3.93±0.14 1.03 

NOTE: (1): Region identification numbers. (2)-(7): Emission line fluxes ratios along with their errors. 

These ratios are measured relative to Ha, where the Ha fluxes are measured in units of 
10~ 16 ergsec -1 A -1 cm~ 2 arcsec -2 and are absolute-flux calibrated. The line fluxes listed are those of: 
[OII]A3727; [OIII]A5007; [NII]A6548; and [SIL1A6717, A6731 (from left to right). (8): Ay is calculated from 
the Balmcr decrement with adopting a MW extinction law (Ay/A^ =1.164, Cardelli et al. (1989)). * The 
term total refers to the integrated spectrum. In this case, the flux-lines ratios values are in units of 

10~ 12 erg sec -1 cm~ 2 . 



Table 4: Ha equivalent widths in emission from the ring spectra. These values are not corrected for underlying 
stellar absorption. 



Ring 


EWgS 
(A) 


Radius 
(") 


Ring 


(A) 


Radius 

(") 


1 


t 


1 


9 


38.34 


40 


2 


8.75 


5 


10 


37.85 


45 


3 


18.65 


10 


11 


36.5 


50 


4 


29.32 


15 


12 


24.04 


55 


5 


25.56 


20 


13 


34.98 


60 


6 


39.18 


25 


14 


50.18 


65 


7 


44.97 


30 


15 


100.1 


70 


8 


44.00 


35 


16 


128.5 


75* 



t No line emission was detected. 

■Tn this table we present the value for ring 16 that it is not shown in any of the other tables in the paper 
since only the Ha line could be detected in the spectrum and its flux (and EW) measured. 
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Table 5: UV-continuum attenuation results. 



Radius 


A-FUV 


A-nuv 


E(B-V)cardelli 


Av,Cardelli 


E(B- V)calzetti 


A-V,Calzetti 


(") 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


6 


1.98 


1.40 


0.18 


0.54 


0.17 


0.69 


12 


1.51 


1.13 


0.14 


0.44 


0.14 


0.56 


18 


1.33 


1.03 


0.13 


0.40 


0.12 


0.51 


24 


1.20 


0.95 


0.12 


0.37 


0.12 


0.47 


30 


1.04 


0.87 


0.11 


0.34 


0.11 


0.43 


36 


0.90 


0.80 


0.10 


0.31 


0.10 


0.39 


42 


1.01 


0.86 


0.11 


0.33 


0.10 


0.42 


48 


1.04 


0.87 


0.11 


0.34 


0.11 


0.43 


54 


1.23 


0.97 


0.12 


0.38 


0.12 


0.48 


60 


1.20 


0.96 


0.12 


0.37 


0.12 


0.47 


66 


0.95 


0.83 


0.10 


0.32 


0.10 


0.41 


72 


0.84 


0.77 


0.10 


0.30 


0.09 


0.38 


78 


0.84 


0.77 


0.10 


0.30 


0.09 


0.38 


84 


0.93 


0.82 


0.10 


0.32 


0.10 


0.40 


90 


0.75 


0.73 


0.09 


0.28 


0.09 


0.36 


96 


0.40 


0.59 


0.07 


0.23 


0.07 


0.29 


102 


0.62 


0.67 


0.08 


0.26 


0.08 


0.33 


108 


0.58 


0.66 


0.08 


0.26 


0.08 


0.32 


114 


0.90 


0.80 


0.10 


0.31 


0.10 


0.40 



(1): Radius in arcsec. (2): Far-UV attenuation values in magnitudes. (3): Near-UV attenuation values in 
magnitudes. (4)-(5): Galactic color excess from UV data and attenuation values calculated via Cardclli law, 
(E(B-V) = Anuv/S-0) Rv — 3.1, expressed in magnitudes. (6)-(7): Galactic color excess from UV data and 
attenuation values calculated via Calzetti law, (E(B-V) = Anuv /8.22) Rv = 4.05, expressed in magnitudes. 
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Table 6: Abundance gradient fits. 



Fits Type 


TBREAK 


Segment 


12 + log(O/i0 at r=0 


(dex/") 


Gradient 
(dcx/kpc) 


(dex/R 25 ) 


Single Slope 

Double Slope Unweighted 
Double Slope Weighted 


40" 
36" 


Inner 
Outer 
Inner 
Outer 


8.640 ± 0.047 
8.876 ± 0.056 
8.269 ± 0.012 
8.949 ± 0.053 
8.269 ± 0.012 


-0.0042 ± 0.0009 
-0.013 ± 0.002 
0.002 ± 0.002 
-0.017 ± 0.002 
0.002 ± 0.002 


0.0346 ± 0.0074 
0.1073 ± 0.0165 
0.002 ± 0.019 
0.140 ± 0.016 
0.002 ± 0.019 


0.416 ± 0.089 

1.29 ± 0.20 
0.198 ± 0.198 
1.683 ± 0.198 
0.198 ± 0.198 



Table 7 

ionized-gas diagnostics in the h ii regions and concetric annuli. 



ID 


12 + log{0/H) 


q 




T e 


Offsets 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Rings 


1 






1137 


53359 


(-3.6,2.5) 


2 






206 


11036 




3 


8.69±0.11 


8.9E+06 


122 


9822 




4 


8.60±0.16 


9.5E+06 


186 


9685 




5 


8.63±0.32 


2.2E+07 


138 


8362 




6 


8.61±0.18 


1.1E+07 


169 


9721 




7 


8.50±0.13 


1.5E+07 


154 


10184 




8 


8.29±0.18 


1.7E+07 


106 


9215 




9 


8.50±0.12 


1.6E+07 


84 


11308 




10 


8.21±0.63 


1.7E+07 


27 


11204 




11 






256 


10124 




12 






248 


15785 




13 






287 


16225 




14 






242 


11186 




15 






56 


10296 




HII Regions 


1 


7.96±0.74 


2.5E+07 


123 


9549 


(-71.9,17.1) 


2 


7.43±1.08 


2.0E+07 




9229 


(-62.8,20.3) 


3 


8.34±0.12 


1.7E+07 




9984 


(-31.5,35.8) 


4 


8.36±0.07 


1.1E+07 




9247 


(-52.1,25.9) 


5 


8.22±0.11 


1.6E+07 




9211 


(-51.7,20.1) 


6 






272 




(-50.7,53.6) 


7 


7.00±1.20 






10372 


(-45.7,38.5) 


8 


8.30±0.21 




405 


9229 


(-40.4,-0.5) 


9 


8.34±0.10 


1.7E+07 




9935 


(-39.9,-7.1) 


10 






10 


10084 


(-38.2,18.9) 


11 


8.38±0.11 


1.7E+07 


68 


11163 


(-34.9,32.2) 


12 


8.33±0.11 


1.8E+07 


33 


11368 


(-34.5,38.2) 


13 


8.43±0.38 


1.3E+07 




11047 


(-31.8,-3.9) 


14 


8.60±0.22 


1.5E+07 




9940 


(-31.3,8.4) 


15 


8.39±0.70 


1.1E+08 


352 


7959 


(-30.9,-28.0) 


16 


8.17±0.08 


2.0E+07 


30 


9175 


(-26.2,29.5) 


17 


8.45±0.36 


1.8E+07 




15032 


(-23.6,-36.8) 


18 








8371 


(-21.6,51.2) 


19 


8.58±0.20 


1.3E+07 


11 


9295 


(-21.1,17.6) 


20 


8.55±0.97 


7.4E+07 


12 


7556 


(-21.0,-69.4) 


21 


8.62±0.09 


1.3E+07 


116 


10033 


(-20.7,8.8) 
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Table 7 — Continued 



ID 


12 + log(0/H) 




q 




T e 


Offsets 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


22 


7.80±1.00 


1 


.1E+07 


142 


10992 


(-18.6,36.5) 


23 


8.18±0.08 


3 


.9E+07 


78 


10008 


(-15.9,-42.8) 


24 


8.48±0.12 


1 


.1E+07 


12 


9657 


(-14.8,47.1) 


25 


8.50±0.23 


2 


.8E+07 


37 


10573 


(-14.2,-58.0) 


26 


8.72±0.25 


1 


.OE+07 


42 


8739 


(-13.6,2.8) 


27 










19479 


(-12.5,77.3) 


28 


8.70±0.15 


1 


.OE+07 


103 


9507 


(-8.9,-6.5) 


29 


8.61±0.26 


1 


.4E+07 


29 


8659 


(-9.2,21.9) 


30 


8.58±0.25 


1 


.4E+07 


194 


8860 


(-7.5,-15.0) 


31 


8.25±0.15 


1 


.7E+07 




9175 


(-5.7,-36.6) 


32 


8.50±0.23 








8699 


(-3.7,-27.9) 


33 


8.69±0.09 


1 


.2E+07 


73 


8941 


(-2.3,-12.2) 


34 


8.39±0.21 


1 


.5E+07 




11198 


(0.8,55.9) 


35 


8.66±0.07 


1 


.2E+07 




7534 


(1.3,23.3) 


36 


8.41±0.04 


1 


.6E+07 




11831 


(1.3,77.6) 


37 


8.48±0.24 


2 


4E+07 


116 


8804 


(1.4,-24.9) 


38 


8.64±0.17 


1 


.OE+07 


50 


9067 


(3.3,-15.6) 


39 










13102 


(3.3,-69.3) 


40 


8.47±0.18 


1 


.3E+07 




13174 


(3.8,46.5) 


41 








168 


8506 


(5.3,-31.0) 


42 


8.57±0.13 


1 


.1E+07 




10684 


(4.8,35.9) 


43 


8.61±0.39 


1 


.2E+07 


54 


7507 


(6.6,20.5) 


44 


8.74±0.05 


9 


.1E+06 


242 


9523 


(7.9,-0.7) 


45 








149 


11738 


(7.8,8.9) 


46 


8.57±0.26 


9 


.1E+06 


126 


8596 


(11.4,-18.7) 


47 


8.51±0.24 


1 


.1E+07 


241 


8870 


(11.5,-25.1) 


48 


8.63±0.12 


8 


.5E+06 


130 


13086 


(13.6,-9.6) 


49 


8.42±0.11 


2 


.8E+07 


170 


10808 


(14.8,68.9) 


50 


8.54±0.30 


1 


.4E+07 


201 


8500 


(19.0,-1.1) 


51 


8.47±0.14 


1 


.7E+07 


118 


9599 


(20.3,25.5) 


52 


8.53±0.27 


1 


.9E+07 


137 


8655 


(23.8,15.9) 


53 


8.31±0.10 


1 


.5E+07 


158 


9683 


(25.0,21.6) 


54 


8.43±0.08 


2 


.6E+07 




10588 


(25.6,59.5) 


55 


8.44±0.08 


2 


.4E+07 


54 


10645 


(30.7,50.5) 


56 








50 


6785 


(33.3,8.9) 


57 


8.57±0.16 


8 


.1E+06 


259 


23652 


(34.3,-0.0) 


58 


8.34±0.21 


1 


.9E+07 


170 


11032 


(42.5,29.0) 


59 


8.34±0.37 


1 


.OE+07 


1318 


12890 


(45.2,13.2) 


60 


7.96±0.83 


1 


.1E+07 


335 


10993 


(48.3,2.9) 


61 








619 


19051 


(69.1,14.2) 
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Table 7 — Continued 



ID 


12 + loq{0/H) 


q N e 




Offsets 


(1) 


(2) 


(3) (4) 


(5) 


(6) 


62 




60 


12023 


(-3.2,4.3) 



NOTE: (1): HII regions and annuli identification. (2): Oxygen abundances along with the corresponding 
errors. (3): Ionization parameter q defined as q = Sho/tt-. (4): Electron density in units of cm~ 3 . (5) 
Electron temperature in Kelvin. (6): Offset coordinates of the HII regions in units of arcsec. 
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